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ABSTRACT 


Condenser is one of the most important components 
in all thermal power plants, large refrigeration and air- 
conditioning units and in many of the chemical and process 
industries. Hence, the use of this equipment embraces quite 
a large arena of the indxas trial systems. The cost of liie 
condenser, therefore, forms a vital part of the total cost of 
the above systems. Any effort to bring about a reduction in 
this cost is an attractive and useful proposition. This thesis 
is devoted to a study in this direction. 

Most of the condensers for commercial purposes, use 
circular tubes for carrying cooling water through them, which 
catises condensation over the tubes. If the tube-profiles are 
changed in a manner, such that, they are drawn in the direction 
of gravity, the effective gravity force, helping in draining 
out the condensate from the tube surfaces will be larger than 
that of the corresponding circular profile tubes . Moreover, by 
suitably adjusting the radii of curvature, so that they 
increase continuously from the upper generatrix to the lower 
generatrix, a favourable pressure gradient is made to develop, 
which also helps in a faster removal of the condensate from 
the surface of the tubes. Hence, a thinner film of condensate 
is left over the tubes having less thermal resistance and as a 



consequence, the heat transfer rate is increased. For a given 
condenser heat-duty, therefore, less surface area is required 
and, thus, there is a reduction in the first cost of the 
equipment. Such tube profiles are named as ' progressively 
increasing radius of curvatiore* (P.I.R.C.) profiles. 

This concept has been explored in the present work. 

t * -y * Q ^ 

Two tube -profiles, given by the polar cu3rves r=a0 and r-as , 
reflected symmetrically about their vertical axes, have been 
used in this context. The effect of using these profiles on 
the augmentation of heat transfer in the condenser in 
comparison with the corresponding circular profiles is also 
studied . 

The first cost of the condenser, however, is not 
representative of the complete picture of its cost. The operating 
or running cost of the equipment should also be considered in 
addition to its first cost. The operating cost of a condenser 
is due to the energy, consumed in the pump for circulating 
cooling water through the tubes. A complete cost analysis is 
also carried out in this thesis to obtain the optimum condenser 
tube dimensions on the basis of minimiam total cost of the 
equipment. A real-life problem is taken by considering the 
condensing unit of the 110 MW unit of 'The Panki Thermal Power 
Station, Kanpur', in order to perform this cost analysis. The 
reduction in size of this condenser and the consequent decrease 



xvi 


in its first cost have also been determined, when circular 
profile tubes of this condenser are replaced by the already 
mentioned P.I.R.C. profile -tubes . 

Since, the solution of the above problem becomes 
intractable by pure analytical methods, numerical techniques 
with the help of a digital computer (DEC-1090) have been 
employed for the solution. Well-doc\amented computer programmes 
in FORTRAN language have been developed for this pixrpose. 
Necessary graphs and the shapes of the P.I.R.C. profiles 
for the condenser tubes are drawn by using the subroutines 
from'lGL' (Interactive Graphics Library). 
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1 .1 INTRODUCTION : 

Condensation occurs when the temperature of a 
vapour is reduced belov.' its saturation temperature. In 
industrial equipments, the process commonly results from 
contact between the vapour and a cool surface . The latent 
energy of the vapour is released, heat is transferred to 
the surface and the condensate is formed. Another common 
mode of condensation is ^homogeneous condensation % where 
vapour condenses out as droplets, suspended in a gas phase 
to form fog. 'Direct condensation' occurs v/hen vapour is 
brought in contact with a cold liquid. Among the modes of 
condensation, surface condensation is the one, which is most 
frequently encountered and hence forms the matter of most 
concern for discussion and study. 

Out of the two different categories of SURFACE 

f t 

CONDENSATION, the one, called dr opwise -condensation has 
aroused more curiosity in the scientific community. In this 
type of condensation, the vapour impinges on the cool wall, 
decreasing its energy and thereby liquefying and forming 
drops of condensate . These drops grow by direct condensation 
of vapour on them and by coalescence with neighbouring 
droplets, until they are swept off the surface by the action’ 
of gravity or other body forces, surface-tension, and/or shear 
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stresses due to vapoui' f],ow. As the drops niove, they 
coalesce v/ith other droplets in their path, sweeping a 
portion of the surface clean so that condensation can begin 

i ' 

anew. The details of dropv/ise condensation are not 
completely understood, but it is known to take place under 
circumstances where the liquid condensate does not wet the 
surface. 'The reported coefficients of heat transfer in 
^dropwise condensation are found to be an order of magnitude 
(10 times) higher than 'filmwise condensation, to be 
discussed later in this chapter. To achieve 'dropwise 
condensation* on condensing surfaces, the surfaces have to 
be coated. Silicones, cupric oleate, several organic polymers 
and an assortment of waxes and fatty acids are often used for 
this purpose. However, such coatings gradually lose their 
effectiveness due to oxidation, fouling or outright removal and 
film condensation eventually occurs . Permanent coatings of 
noble metals (gold, silver, rhodium, palladium, platinum) 
or teflon overcome this difficulty. The substantial increase 
in cost, hov/ever, far outweighs the advantages, gained there 
from. Although, 'dropv/ise condensation' of steam has been 
studied for over fifty years, this unique area of heat 
transfer has remained as a laboratory curiosity, primarily 
because permanent hydrophobic coatings have not been developed 
to the satisfaction of condenser designers . In a very recent 
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paper £^ 2 . 0 ^ 7 , it has been presented, that organic coatings 
are successful in promoting good quality tiropv/ise condensation' 
for periods upto around 12 000 hours. For a continuously 
running pov/erplant condenser, this requires a change oi 
coating, rougrily, after every 1,5 years. This is an 
attractive proposition. But the uncertainity in this prediction 
by the authors has been reflected in their pronouncement, 

*'’The possible deterioration of the heat transfer with 
exposure time to steam remains to be determined. So this 
sends us back to the original impasse. If, some day, the 
technique of ' dropwise -condensation ' is achieved to the level 
of perfection, the dream of ' steam cars ' will become a reality . 
At present, however, it seems to be a little distant reality . 

Since, it has been difficult to sustain ' dropwise 
condensation' commercially for long periods of time, conser- 
vative designs are based on operation with the other signi- 
ficant mode of ' surface condensation,' called the filrawise 
condensation.' In this mode of condensation, the liquid 
condensate forms a continuous film, which covers the surface 
and takes place when the liquid wets the surface . This film 
flows over the surface under the action of gravity or other 
body forces, surface-tension and/or shear stresses due to 
vapour flow. Heat transfer to the solid surface takes place 
through the film, which forms the greatest part of thermal 


resistance . 
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Condensation heat transfer can be studied from 
either the macroscopic or microscopic point of view, as is 
the case in any other physical phenomenon. 'The ultimate 
objective of the engineer is to understand the phenomenon 
so that it can be described precisely and hence, its behaviour 
can be predicted for design purposes. The macroscopic view 
point considers matter as continuum and requires the intro- 
duction of certain phenomenological laws, the macroscopic 
transport equations, in addition to the conservation equations 
of mass, momentum and energy. On a microscopic level, dealing 
v/ith molecules, atoms and subatomic particles, only the 
conservation principles need be considered, Hov/ever, because 
of limitations in describing these particles and in relating 
the large nimiber of particles, which constitute a reasonable 
size system, one resorts to the concept of statistical mechanics. 
For practical applications to engineering systems, however, 
only the macroscopic method of study is ■ of relevance. 

Most of the extensive research in the field of 
condensation heat transfer has been devoted to the understanding 
of the basic process. In recent years, however, there has 
been an increased emphasis on techniques to augment this two- 
phase heat transfer process . Efforts are continuing to devise 
means for predicting the heat transfer coefficient for 
condensation more accurately and to increase the heat transfer 
coefficient. This will result in reduced condenser sizes and, 
hence, the first cost of the unit will be reduced. 
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The increase in heat transfer coefficient will result, 
if ' dropv/ise condensation' can be sustained, or if the thickness 
of liQ_uid film, in case of ^ film condensation, can be reduced. 
The first of these alternatives is of no practical interest, 
at least for the present, in the light of ^^/hat has been already 
discussed on the topic of 'dropwise condensation 'in this 
section. Hence, almost all practical methods, used to augment 
condensation heat transfer, aim at achieving the second 
alternative i.e. reducing the thickness of the condensate 
film on the condensing surface. In most of the industrial 
applications, e,g. large thermal power plants and refrigeration 
units, condensation takes place in a closed shell over tubes, 
that carry the coolant through them. Fig. ( 1 . 1 ) shows the 
outline of such an arrangement, as used at the PANKI THERMAL 
POV/ER STATION, KilNPUR.' Such a process of condensation is 
termed in technical jargon as 'vapour space condensation' or 
'external condensation.' 'The application of condensation inside 
tubes, called ' In -tube condensation,' enjoys very limited 
applications like domestic refrigerators and is, therefore, 
not discussed further. 

The above discussions lead one to conclude that from 
the engineering viewpoint, the study of film condensation by 
macroscopic modelling, forms the most important basis for 
studying the phenomenon of condensation. 
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^ ‘2 flEVIEW OF PREVIOUS ;VORK : 

Most of the condensers, used in commercial practice 
are designed on the basis of ^filmwise condensation/ The liquid 
film, formed over the condensing surface, forms a significant 
source of thermal resistance in the path of heat transfer from 
the vapour -to -coolant heat exchange process. Any attempt to 
reduce the thickness of this liquid film by augmenting the 
condensate drain-out rate, formed over the condenser tubes will 
therefore, reduce this thermal resistance and, v/ill, thus, 
enhance the heat transfer coefficient of the vapour -side. 

Much work has been done in this field, dealing with 
the augmentation of film condensation process and many 
research papers have been published over the years. A.E. Bergle 
_^14,15_7 has made a very useful review of most of this 
ever-widening literatux'e. lie has deal L with these techniques 
of augmenting condensation heat transfer' under the categories 
of (i) active techniques and (ii) passive techniques. The 
active methods are those which I'equire external power to 
effectuate the process such as the use of mechanical aids, 
surface vibration, fluid vibration, electrostatic fields, 
^injection and suction etc. Passive methods, on the otherhand, 
require no such external power. Some of the passive techniques, 
used successfully include treated surfaces, rough surfaces, 
extended surfaces, displaced enhancement, swirl -flow devices, 



surface-tension devices and liquid additives. From the 
point of view of what has been done in the present thesis, 
the development of surface-tension devices as a passive 
technique, needs a little elaboration. 

• Gregorig /~1o,l6_7 was the first to point out the 
possible advantages of using surface-tension-driven-forces 
in enhancing the c jnuen.sin,; heat trin;; I'er I'.oo f i'icien t . His 
concept of transforming Uiis idea into practice materialiiied 
in the development of what is Imovai as the [’luted tube, as 
shovm in Fig. ( 1 . 2 ). This is essentially a tube with fine 
corrugations on the surface. The curvature of the surface of 
the condensate film, formed on a fluted tube is such that 
the condensate is driven by surface -tens ion -induced forces 
from the ridge to the trough. The mechanism, responsible here 
is that the curvature of the film-surface on the ridge side 
gives an increased pressure there, while the pressure on the 
trough side is reduced by reversing the curvature. Hence, 
there is a pressure gradient from the ridge to the trough, 
which helps in driving the condensate from the ridge. This 
leaves a very thin film over the ridge with low thermal 
resistance. Thisassists in the augmentation of heat transfer. 
Since, the 'fluted tube' is very complex in shape, incorpora- 
tion of this device in heat exchanger tubes would require 
a detailed a-priori economic analysis. 




FIG. 1.2 PROFILE OF THE CONDENSING SURFACE 
DEVELOPED BY GREGORIG . 



A relatively simple concept for augmenting 
condensation heat transfer has been pronounced by G.G.Shklover 
et. al. ri_7. In this paper, the authors have conceptualised 
the use of tv/o very important factors, v/hich help to 
enhance the condensate drain-out rate from the condensing 
surface. Circular profile tubes are moulded into the shape, 
as shovm in Fig. (2.1.b), keeping the perimeters constant. 
Since, the tubes are drawn more in the direction of gravity 
in comparison to the circular ones, the contribution of 
gravity force will increase the drain -out rate of the 
condensate. Moreover, if in the changed tube profile, the 
radius of curvature is made to increase continuously from 
the upper generatrix to the lov/er generatrix, a favourable 
surface-tension pressure gradient develops, which assists 
in the further enhancement of the liquid drain-out rate. 

Both these effects viz, (a) Gravity force and 
(b) surface-tension force, in unison, help in the augmenta- 
tion of condensation heat transfer. In spite of the concept 
being excellent, the paper does not incorporate a clear 
and complete solution oT the problem, llov/evcr, a suggestion 
has been given by the authors that apialica tion of this 
concept to real-life problems, such as condensers in steam- 
power -plants , will be quite useful. The profile of the tubes, 
chosen by the authors is that of a logarithmic spiral, 

v/hich is a little complex, the pol/ir form of which is given 
Ocota , 
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No reviev/ of li ter-i ture in the aren oi film- 
condensation v/ill be complete v/ithout mentioning the pioneering 
paper of Nlisselt v/hich was pui.lished in the year 

19'16 and still remains the only original paper to describe 
the phenomenon of film -condensation analytically. 

Improvements have been made over the years, but excepting 
the behaviour of liquid metals, the original tl^ory has been 
reasonably successful, ilven the authors in have used 

the same analytical approach as was used by Nusselt 
to formulate their problem. 

1,3 PRiSSENT WORK : 

The present work is mostly an extension of the 
paper £~ 1 _ 7 . A detailed and complete solution is obtained 
for the problem. The suggestion of the authors in /~ 1_7 is 
also incorporated. The complete analysis is made on the 
basis of data, collected from the " PANKI THERML POV.'ER 
STATION, KANPUR for the condenser of their 110 M//. unit. 

Two tube profiles, represented by the polar curves 
^ r=a9' andV=ae®^, respectively, are considered to take 
advantage of the augmenting effects, mentioned in Section 
(1.2). The economics of this change is also discussed. 

Employment of the above mentioned tube profiles, 
however, is against the conventional practice of using 
circular tubes in condensers, v.'henever, a change is 
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contemplated in the conventional practice, more often than not, 
there is a tendency to view it with an attitude of inertia . 
Hence, the use of changed tube profiles in the condenser 
may be considered as a long-term solution towards a 
rational designing of a compact heat emchan, er. 

A short-term, useful approach Ls also adopted, 
not from the point of viev/ of reducing the condenser size, 
but on the basis of a cos t-optimiscation technique. A brief 
description of the concept is as follov/s : 

A reduction in the tube diameter's of the condenser 
results in augmentation in the overall heat transfer 
coefficient. Hence, for a given heat transfer requirement 
of the unit, less tube material is required and the first or 
procurement cost of the unit is reduced. But in consequence, 
the pressure drop through the tubes is increased and this 
requires a higher pumping power. This entails a higher 
pumping energy cost. Both these effects work in opposite 
directions, as far as the objective of minimizing the 
total cost of operation is concerned. Keeping in view these 
two contradicting requirements, the optimum tube -dimensions 
are determined on the basis of minimum total cost. 


The same analysis is also extended to the 
changed tube profiles with suitable modifications . 
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Computer pror^rammes in FOUTHjIM langua;:;e are 
developed for the solution of all the above problems . 

The programmes are made very general to encompass a 
variety of similar problems, ’.vhich can be solved v/ith 
very little modification in the parent programmes. They 
are properly documented for easy understandability of any 
user, having a working familiarity with FORTFlAil language. Mo 
subroutines (such as those, supplied by 'Numerical Algorithm 
Group/ abbreviated as ) are used in the programmes. 

Hence, in those cases, where accessibility to sophisticated 
computing systems like DEC-10 is not available, IBM personal 
computers with FORTFbiM facility can be used to obtain the 
solutions v/ith the help of the above programmes . Interactive 
Graphics Library/ popularly known by its acronym IGL 
subroutines are used to draw graphs and the shapes of the 
tube profiles. 
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--*I'ii\LYSI3 OF COiiDENoATIOM HEAT T 
tation by using Pit ogress IVELY I 

OF CURVATURE PROFILE TUBES . 


”■<' -R ‘P / ' J ■' T T r'*- rj-T*.'- 

:rii.ASING RADIUS 


2.1 INTRODUCTION : 


One of the major heat transfer resistances, 
associated with the process of steam, condensing over 
horizontal tubes of a condenser, is that of the condensate 
layer, formed over the tubes. Hence, heat transfer with 
film condensation is limited by the thicimess of the 
condensate layer. She growth and removal of this condensate 
layer ciepends on the forces, acting on the film. For a 
tube v/ith circular profile, which is the most general in 
applications, the velocity of the condensate is solely 
determined by the projection of the gravity force along the 
tangent to the profile and the viscous force. 


If the profile of the above circular tube is moulded 
in such a manner that it is drawn in the direction of the 
gravity vector (Fig. 2.1(b)), the average value of the 
projection of the gravity force along the tube periphery 
will be greater than that of the corresponding value for 
circular tube . This increases the rate of run-off of the 
condensate and the heat transfer coefficient is, 
consequently, increased. 
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In addition to the above modification, if the 
moulded tube profile is made to have its radius of curvature, 
continuously increasing from the upper generatrix to the 
lower generatrix, then, since the surface-tension force is 
inversely proportional to the radius of curvature, a 
pressTire gradient develops along the tube surface. Since, 
the pressure constantly decreases towards the lower 
generatrix, it acts favourably to enhance the condensate 
drain-out rate. The circular profile, having thie same radius 
of curvatxire throughout, does not develop any pressure 
gradient and, thus, cannot take advantage of the effect of 
a favourable pressure gradient. 

Hence, this discussion reveals the fact that the 
P.I.R.C, (progressively increasing radius of curvatxire) 
tubes (Fig. 2.1 .(b)), by making use of the increased 
gravity force and a favourable pressure gradient (due to 
surface-tension), result in a faster condensate drain-out 
rate from over the tube surfaces, when compared with their 
circular counterpart. Consequently, the heat transfer 
coefficient registers an increase and for a given condenser 
heat transfer requirement, a relatively smaller heat 
exchanger can serve the purpose. 
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In this chapter, details of the solution method 
are presented for three different tube profiles j viz.; 

(a) Circular 

(b) Polar form V=aO', being reflected on both 
sides of its vertical axis . 

{Cj Polar form r=ae , being reflected on both 
sides of its vertical axis. 

The solution for circular profile is already 
available from Nusselt’s analysis. 'The present solution, 
then, is compared with the available result to vindicate the 
impeccability of the method, employed. Profiles (b) and (c) 
are chosen in line with the above discussion to take advan- 
tage of the gravity and surface-tension effects to augment 
the heat transfer rate in the condenser. The percentage 
increase in the average steam-side heat transfer coefficient, 
using (b) and (c) in comparison to (a) is determined. The 
isolated and combined effects of gravity and surface-tension 
on the augmentation of condensation heat transfer coefficient 
are presented, separately. The savings in first-cost, 
accruing from a smaller size condenser using I'.I.R.d. profiles 
(b) and (c) in comparison v/ith the presently used circular 
profiles are also determined. 
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^ ^ ) i-lADii IN THil Ail-iLYoI6 'H-:-^iK v .■-.Llijl'ry : 

The actual process of condensation is very 
complex. In order to make the problem of film-condensation 
solvable by analytical technique, certain simplifying 
assumptions are made. These assumptions make the solution 
of the problem tractable and still preserve its realistic 
nature. Despite the complexities associated with film 
condensation, useful results are obtained by making the 
following assumptions . 

1 . The flow of the condensate film is laminar. It has 
been found practically /~18_7 that the turbulent flow of 
the condensate is difficult to achieve on horizontal tubes, 
but may be easily established over the lower portions of 
vertical surfaces, ■-^ince, in this analysis, condensation 
over horizontal tubes is considered, the assumption is 
sufficiently reasonable . 

2. The steam to be condensed is at a uniform tempera cure 
T , . There is no temperature gradient in the vapour 

boundary layer, formed over the liquid condensate layer . But 
in actual practice, a temperature diffei'ence in the 
vapour boundary layer is observed to manifest . However , 
its magnitude is of the order of 0.03*^C(for steam), whicn 
for all practical purposes can be neglected. Moreover, it 
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has oeen observed. that the thickness of the vapour 

boundary layer is of the order of 2 to 3 times the molecular 
diameter. For steam, the thickness is of the order of only 
one molecular diameter. The consideration of these microscopic 
dimensions in an analysis, based on the continuum approach, 
loses significance. Both these effects provide sufficient 
ground to validate the assumption, flence, heat transfer to 
the liquid -vapour interface can occur only by condensation 
at the interface and not by conduction from the vapour. 

3. The shear stress at the liquid -vapour interface is 
assumed to be negligible . '%is assumption is true because of 
the low steam velocities (< 1 7(^m./sec . ) , commonly encountered 
in Vapour space condensation /~21_7 which prevents the 
interfacial drag from becoming of any consequential importance, 

4. T};ie convective terras in the momentum and energy 
transfer in the condensate film are assumed to be negligible . 
This assumption is true because of the very low magnitude of 
velocities ( of the order of 10 m. /Sec.), encountered by the 
condensate film. Such low velocities are inconsequential from 
the point of view of engineering applications. Hence, the 
boundary layer flow of the condensate is one of very slow 
motion, commonly termed as the 'creeping motion.' 
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5. The surface temperature of the tube on the steam- 
side is assumed to remain constant. Since, the magnitude 
of heat-transfer coefficient of condensation on the steam- 
side is very large, the relative order of magnitude of this 
coefficient requires an almost negligible temperature 
difference for heat transfer. Moreover, actual experimental 
observations by Peck and Reddy and by Bromley, Brodkey and 
Fishman separately, corroborate the validity of this 
assumption . 

6. Alongwith these assumptions, Nusselt's modification 
of the analysis for condensation over 'n' vertical tubes in 
a row and Chen’s modification for the sub-cooling of the 
condensate film are also assumed to apply. This assumption is 
valid for those cases where Jakob No. 'Ja^ is less than 2. 

In the present case, 'Ja’ has been calculated out to be 
0.0035 and thus speaks for the validity of the assumption. 
This assumption is made to make the analysis more 


realistic 
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2.3 THE TUBE -WALL SURFACE TEMPErlATURE ; 


The average steam-side heat transfer coefficient 
for circular profile is given by the well-known relation /~5_7, 


E= (0.728)x 


"p-, ('^sat.-^w) 
1+0.2 ir (n-1) 


h 


fg 


*2 3 ' 

PlS ^1 ^fg 




5r^7iirrr:+ -'C) 

o sax. 


( 2 . 1 ) 


The above relation is valid when film-condensation 
of a vapour is taking place over 'xi' rows of vertical tubes 
of a condenser and the effect of condensate sub-cooling is also 
considered. 


In expression (2.1), the properties , 'p^' 
and are determined at the reference temperature 

'^ref. where, 



^ref. = \ " 0*25 (T^^^^-T^) 


( 2 . 2 ) 


Knowledge of the tube-wall surface temperature 'T/ is 

fV 

essential in order to determine the heat transfer coefficient 
'£g. The following procedure is adopted to obtain 'T^' . 

Considering a set of two tubes (one each for entrance 
and exit of cooling water) of 7,5m, length* each, which are 
presently used in idie condenser, taken for the case -study and 
using 'Appendix-1 ' for the operating data of the above mentioned 


* The procedure is applicable to any general ^.ength of the 
condenser tube. 
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condenser, the heat transfer rate ' Q ' is given hy, 

Q = 21890 W 

Also, 'Q' is expressed as. 

In the above expression, 'A^ and ' can be obtained 

using Appendix -1 . 

Expressions (2.1) to (2.3) for ' "^ref '» 

all depend on It is, therefore, necessary to first 

determine ' T ' . This has been achieved as follows : 
w 

Variation of the thermiophysical properties 'P-i 'C ' , 

X 

' k£ and ' jul^' with temperature are tabulated and given in 
'Appendix -2 A 'computer programme No.1 ', which uses the 
method of least squares , has been developed to obtain 
polynomial relations for these themiophysical properties in 
terms of the condensate temperature. 'Appendix -3' lists the 
/computer programme No.l ' and ' Appendix -4 ' gives the 
polynomial relations. These polynomial relations have been 
substituted for the thermophysical properties ' in expression 
(2.1) for determining 'S' . A 'computer programme No. 2', 
listed in 'Appendix -5' has been developed as an iterative tool 
to evaluate ' T^'. 



Initially, a guess value of ' is fed into the 

'computer programme No,1 ' . With this guess value, 'h' is 

s 

determined using expression (2.1.). '0' is then calculated 
by using expression (2.3). This value of ' Q' is then matched 
with the actual magnitude of 'Q', which is 21890 W. The 
difference between the two is noted. A new guess-value is 
assigned to 'T^ and the procedure is repeated till the 
difference between the actual and calculated values of 'Q' 
is negligible. 

'Table 2.1' gives resxilts for 'T^' and some of 

the related properties of the condensate, corresponding to 

this value of ' T ' . 

w 



2.^ general formulation of the problem to DETERI4II^ 
CONDENSa^IOIS^ HSA^ T^^tlSI'ER C^EF^lClJ^NT ; 


The problem of determining the heat -transfer 
coefficients for the circular and polar-profile -tubes requires 
specifying the governing differential equations and the 
necessary boundary and/or initial conditions. The coordinate 
system used is shown in 'figure 2,1'. The effort, here, is to 
first derive expression of the average heat transfer coefficient 
for the circular tube by solving the general problem and 
comparing the results, so obtained, with the well-known 
correlation (2,1) for 'S'. The same analysis will be extended 
to derive expressions for average heat transfer coefficients 
for the polar-profile-tubes. 

The general momentum and energy conservation 
equations, in their differential form, are, respectively, 
given as following ; 

^1 = +X (2.4) 

h (=p t) + I- I (Cp^T )_7= - p( l-V + I 

+ 0 + Qg (2.5) 

Using the assumptions of steady flow, creeping 
motion and Prandtl's laminar boimdary layer approximations, 
the momentum equation reduces to , 



2 ^ 



3 u 




X 


( 2 . 6 ) 


The energy transfer across the condensate -film 
occurs only by conduction as discussed in Section (2.2). 
Moreover, with the assumptions of 'steady' and 'creeping ' 
motion for the condensate-film, expression (2.5) boils 
down to , 


0 = V . (k^ V T) 


(2.7) 


Assuming constant property of the condensate and using 
Prandtl^s laminar boundary layer approximations, the final 
form of energy equation is given by, 


X 



0 


where, is the conductivity of the condensate and is assumed 
constant . 


3^T 


( 2 , 8 ) 


Solution of the momentum equation for obtaining the velocity 
distribution in the condensate boundary -layer, requires the 
following two boundai*y conditions ; 
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u(y) =0, y= 0 (Condition of no slip ) 

3 

3 y u(y) = 0, y=6 (Zero interfacial shear stress ) 

Similarly, for the determination of temperature distribution 
in the condensate boundary layer, the energy equation (2.8) 
is used alongwith the following two boundary conditions : 

T = T^, y=0 (Constant wall -surface temperature) 

T = '^sat.* (Constant vapour-liquid interface- 

temperatxjre ). 


Since ' u' and 'T' are functions of 'y' only and 'p' 
is a function of 'x' only, the partial differentials are 
replaced by total differentials and the formulation of the 
problem is given by : 


(1) The momentum equation ; 


^ - 1- se + n Y 

dy2 - 


0 


(2.9) 


Boundary conditions : 

u = 0, y = 0 (2,10) 

* 0 , y=6 (2.11) 

(2) The energy equation : 
d^T 



0 


' ( 2 . 12 ) 



Boundary conditions : 
T = \. y = 0 


(2.13) 


=Tsat.’y=^ 


(2.14) 


Expressions (2.9) to (2.14) represent the complete 
formulation of the problem. 


2.5 SOLUTIONS FOR CONDENSATION HEAT TRANSFER COEFFICIENT 
FOR SPECIFIC TUBE -PROFILES ! 

In the absence of any specific method known for the 
solution of the above mentioned problem, a numerical 
technique has been developed for the solution of the problem 
and is tested for the circular profile case for which the 
results are known. The method is then, extended to determine 
heat transfer coefficients for the polar profile tubes. 

2,5,1 Circular Profile : 

The velocity distribution in the condensate boundary 
layer is obtained by solving the momentum equation. In the 
generalised momentum equation (2.9), the following modifica- 
tions are made . 

(a) Since, the surface-tension contribution is absent for 
circular profiles, its presence in the momentum 
equation is dropped, 

(b) '%e body -force component is contributed by gravity only. 
Resolving the gravity vector along the tangent to the tube- 
profile (Fig.2.1), the magnitude of the body -force ’ 

is given by. 
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^ = P]_g Sin P 

= g Sin(0 -90°) 
= Pi g (-CosW; 


^ = P^ g (-CosQ) (2.15) 

v/ith the above modifications, the expression (2.9) becomes 

— ^ (-Cos 0) = 0 (2.16) 

dy*^ 



fi£ 


(-Cos 9) = F(0) 


(2.17) 


Then the momentum equation becomes 

^ + F(9) = 0 
dy^ 


( 2 . 18 ) 


Integrating the equation (2,16) twice and using the 
bo^Indary conditions, given by expressions (2.10) and (2.11) 
the velocity distribution is obtained as : 


u 


Ppg 


(-Cos 9) 


/-a y 


2 

h 


7 


( 2 . 19 ) 


The mean velocity ' u^ at any location along the periphery 
of the tube is then given by, 

J u dy 
o 


1 

= 3- 


( 2 . 20 ) 
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Where, '6' is the condensate -boundary layer thickness. 
Simplification of expression (2,20) yields, 


u 


m 


P]_g 




( — C os 


&) 


( 2 . 21 ) 


The temperature distribution is, then, obtained by 
integrating equation (2.12) twice and using the boundary 
conditions, given by expressions (2.13) and (2.14). This gives 

T = - .y + T^ (2.22) 

Consider now a control volume of imit thickness, 
indicated by the hatched portion of Fig. (2.1 .a), which is 
contained by the infinitesimal angle M0 

An energy balance of this control volume gives : 
dm X hfg - ^ X (dx X 1) (2.23) 

Where , 

dm = the net mass efflux from the control volume. 

= d /~Pl % J (2.24) 

Substitution of the expression (2.21) for in (2.24) 

yields, 

2 -3 

Pt g ^ 

dm = d - I— (-Cos ©)_7 (2.25) 



?0 


dx 


= Latent heat of evaporation 

= Elemental peripheral length 
= R d9 , ’R' being the (2.26) 

outer radius of the circular profile tube. 


But, 



(2.27) 


v/here, 'd^* is the outer diameter of the circular tube. 
Using (2,27), the expression (2.26) for 'dx' reduces to : 


dx = 



(2.28) 


It is to be noted that in expression (2,25), the 
boundary layer thickness ^ d' is a function of ' © 


Substituting for 'dm' and 'dx' in (2,23), we get. 


Pi g 

r (-Cos 0)_7x hjg 


(2.29) 


Simplification of (2.29) gives. 


[T ' I c — TTT 

sat. ^ 


1 


d© 


(2.30) 


In the above expression, the condensate properties are 
evaluated at reference temperature ' 


by using Aptjendix- 4 
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If, 



and 

(- I Cos 9) = F^(9) 

the equation (2.30), may be written as, 


(2.31) 


(2.32) 


(N) X d F^(0). 6- 


. d© 


(N)x 6 X d F^ (9) 6^ = d© 


(2.33) 


Multiplying both side of expression (2.33) by ^F^(©)| 
and simplifying, the following result is obtained : 


^ L 

(N) X d ^ X 6^x F^(©) 


= iFl(9) d© (2.34) 


Integration of equation (2.34) and simplification gives ; 



(2.35) 


The steam-side heat transfer coefficient 'h' at any location 

s 

‘©^ is given by [19], 
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(2.36) 


Substitution of expression (2.35) for '6' in expression 
(2.36) gives ; 



(2.37) 

Using expression (2.31) for 'N' in (2.37), we get, 



(2.38) 


Relation (2.38) is derived, for steam, condensing over 
a single circular, horizontal tube. But in an actual condenser 
condensation occurs over a number of tubes and also the 
condensate is subcooled. To take accoxmt of these effects. 



33 


i'-iusselt’s correction for condensation over ' n' horizontal 
tubes in a vertical row and Chen’s correction factor for 
condensate subcooling /5_7 are applied to expression (2.38). 
The application of these correction factors makes the analysis 
more realistic, '^ith these modifications, expression (2.38) for 
’h^ reduces to , 



III 



It is to be noted that in the expression (2,39) for 'h^ , 
part-I is a function of '9% whereas part-II and part-III 
are constants. Evaluation of the integral in expression (2.39) 
analytically, is a formidable task. In order to overcome this 
difficulty, the following method is employed : 
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The serai -perime ter of the profile is divided into 

100 equal parts. At each discrete location of 'O', the local 

heat transfer coefficient 'hi is evaluated with the help of 

expression (2.39), the integral being obtained by using 

'Simpson^s numerical rule of integration'. Since at '9 =iy '2 ' 

and 0=3%^ the expression (2.39) assumes the indeterminate 

(zero/zero) form, these values are suitably truncated to 

circumvent this difficulty, for example, '9= J 'is truncated to 

~ +0.01jt = 1 .6022' and '9 = ' is truncated to 

^ ■“ = 4.6809' . A 'computer programme No. 3' is 

given in 'A.ppendix -6', which has been used to calculate 'h' at 

s 

discrete locations '9'. Some representative values of 'h' vs.' 9', 

s 

obtained from this computation are presented in ' Table 2.2 '. 

V/ith these values of ' h' and ' 9 ' at discrete locations, a 

O 

polynomial relation between the two variables is obtained by 
using 'the method of least squares ' . A ' computer programme ' 

No, 4 , listed in 'Appendix-7', has been used to obtain the 
polynomial relation between 'h' and '9'. The constant terms, 
viz. part-II and part-III in expression (2.39) are kept unaltered 
and the part-I of this expression (2.39) only generates the 
polynomial relation, which is computed, using the above 
mentioned procedure . 


This polynomial relation between / hi and '9 ' comes 
out to be , 


h„ = 


1 .3324 - 0.4932 © + 0.19^6 


e^-0. 


0293 0 


3‘ 


K 


(2.40) 
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K = 




Pn g K 


^1 


1 

1/4 i 




(2.41) 

Hence, 'K' is the product of the two constant terms, viz. 
part-II and part-III of expression (2.39). 

'Ihe average steam-side heat transfer coefficient ' fi ' Is then 

s 

given by, 

‘'s' Topep— (2.42) 

Substituting expression (2.40) for 'h' in (2.42) and evaluating 
the integral between the truncated limits ' G..j=1 .6022 ' and 
^0^= 4.6809' , we get the desired expression for the average 
steam-side heat transfer coefficient, 

Eg = (0.7303) X [k] (2.43) 

where ' K' is given by expression (2.41). 


The standard result for ' , for circular profile, given by 

expression (2.1) may be represented as ; 

Eg = (0.728) X [k] (2.44) 

Comparison of the derived expression (2.43) for 'E 'with the 

s 

standard result (2.44), it may be easily seen that the percentage 
error between the calculated and known expressions for ' h^t in 
case of circular profile is given by, 



*/. 


error = 


'S' from expression (2.43) -'n' from expression (2,44) 

S * S 


'h ' from expression (2.44) 


-x100 


(0.7303)x [k] - (0.728)x [k] 
(0.728) X K 


X 100 


_ 0.7303 - 0.728 „ 

" 07758 ^ 

= 0.288'/. 


Since, the error is negligible, the calculated result 
is in close agreement v/ith the known result. This proves the 
efficacy of the method, employed for the solution oi the 
problem , 

2,5.2 Profile Represented by polar curve, 'r '» Reflected 
symmetrically on both’~sides of its vertical axis : 

Fig. (2. 2. a) shows the actual shape of the profile 
'r = aO.' The corresponding circular profile, which has been 
moulded to generate the above polar profile is also shown in 
Fig. (2.2.b) . A 'Computer programme No. 5 , listed in Appendix-8, 
has been used to draw these profiles. The same co-ordinate 
system as that of Section l2.5,1) is used. 

Some geometrical relations, concerning the profile are 
evaluated first, which aid in the formulation and solution of 


the problem . 
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SEMI PERI!4ETER 


A differential length 'dx ' along the periphery 
of the profile is given by, 


Here, 




r = a© 
dr 


d© 


= a 


(2.45) 

(2.46) 

(2.47) 


Substitution of (2,46) and (2,47) in (2.45) gives, 

(2,48) 


dx = a ^ / 1 + 0' 


d© 


The semi perimeter 'S^ is represented as : 




= S 
*=”/2 


a /1+ 9 d0 


or 


* i 


9 




On simplification, expression (2,50) gives : 
S = 10.3985 a 

Jr 


(2.49) 
©=3n ! 2 

®=V2 

(2.50) 

(2.51) 


* Results of standard in tegrat ion table have been used 


'x^" H-a'^ dx = ^ X ^/a^+x^ + -^ Sin h”^ (x/a) 


Sin h"‘^(x/a) = In |^x+ -J x^+a^^ 
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Since, this seniip-. rimeter is kei)L equal to Uie 
semi perimeter of the coi'respondln}/ circular profile, il 
leads to the following expression : 


TtR = Sp = 10.3985 a 

^ = 10.'5§85 ^ = 0.3021 R (2.52) 

In the above expression, the factor '0.3021' is 
named as ' radius multiplying factor' or 'RMF '. Multiplying 
' RMF' by the radius of the circular profile, results . in the 
determination of the characteristic constant 'a' of the polar 
profile, ' r = a© 


AREA OF CROSS-SECTION 

The cross-sectional area of the circular profile, which 
is moulded into the polar profile under consideration, is 


given by , 

^■^cn^c 



(2.53) 


where , 

subscript 'c' stands for 'circular.' 

d^ = Inside diameter of the circular tube 

A = Area of cross-section 
cn 

^rom the relation (2,52), the constant a can be 
expressed as : 

d. 

a = (0.3021) X ( ^) 


(2.54) 
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The area of cross -section of the polar profile is given by, 

O D 

= J I- dS 

®=V2 


3tx 


=a 


Iz 


0^ d© 


©=n 


fz 



3n/^ 

®=V2 


= (0.3021)2 X ( X ( ^)2 

3 d 

(0-3021)2 ^ ^ ^j2 (2.55), 

In expression (2.55), the subscript 'p' stands for 'polar.' 

An area factor ' AF ' is now defined, which, when 
multiplied by the cross-sectional area of the circular profile, 
yields the cross-sectional area for the corresponding polar 
profile. In the present case, 'AF' is calculated as : 


AF 


^^cn^pi expression (2.55) 
expression U.^3) 


0.9758 


(2.56) 


HYDRAULIC MEAN DIAMETER 

The perimeters of the circular and polar profiles are 
kept equal to provide the same heat transfer surface area for 
both the profiles. Hence, the perimeter 'PHR' of the polar 
profile, ^ r=a©,' is given by, 
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PER = Ti 


(2.57) 


The hydraulic mean diameter ' HIO ' is then obtained 


using : 


4 X AF X (A ) 

™ m ^ 


(2.58) 


RADIUS OF CURVATURE 


The radius of curvatxxre 'R ' of the profile at any 

c 

location, is evaluated as follows : 




rA 2( ||)2- r(2% ) 


d©' 


(2.59) 


Here, 


„ „ dr „ d^r p, 
r = a9j 3^ “ ^ > — 7 ~ 0 

d©^ 


R 


(a^©^ + 

° a^©^ + 2a^ 


or, 


R„ = a 
c 


(2+©^) 


3/2' 


(2.60) 


SURFACE-TENSION PRESSURE 

The surface-tension pressure at any location is 
inversely proportional to the radius of curvature at that 
location and is given by. 
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P = 


Rq a 


(2 + 0 ^) 


(2.61) 


where , 


= Siirface-tensioii force per unit length. 

The pressure gradient along the profile is, then, 


determined as given below 


( iS) / (i) 


(2.62) 


Using (2.48), 


do = a yi+ 0 


(2.63) 


Differentiating expression (2.61) with respect to ' 0 ^ 
and v/ith the help of (2.63), the expression (2.62) gives ; 


‘"I 

T 


(Q^-t-4Q 
(1 + 0 ^) 


( 2 . 64 ) 


The negative sign of this pressure gradient is in 
conformity with the expectation of a favourable pressure 
gradient . 


Having completed the eval\jation of the above relations, 
we now proceed to the formulation and solution of the problem. 

The governing equations of the general formulation, 
given in Section (2.4) are now modified as following : 

(a) The gravity force is again the same as that for the 
case of circular profile, i.e. ' X - g(-Cos 0) 
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(b) The surface-tension pressure gradient is now given 
by the expression (2.64), i.e. 


^ 1 
3 $ = " 72 

a 


(e^+ 4©) 

X " ' T ^' 

( 1 + <^Y 


Using (a) and (b) in the momentum equation (2.9) » 

we get, 

Momentum equation : 


d^u , 
dy^ ^ 


( —Cos ©) + 



(©^+4©) 
(1+©^^ 


zs 0 


(2.65) 


Boundary conditions ; 


u = 0 , y = 0 



Energy equation is the same as that for the case of 
circular profile. It can be easily seen that in the above 
formulation of the problem for polar profile, 'r=a6', everything 
is the same as the formulation for the circular profile 
(Section 2.5.1) except the momentum equation (2.65). Hence, 
the momentum equation is suitably modified to fit into the 
pattern of the formulation for circular profile in the following 
manner : 

Rearrangement of tenns in expression (2.65) gives : 
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d^u 


(Cos 0) 



(0^+40) ( 

(1+©^)^ j 


consider, 

Fg(e) 


Pi s 

-rr (-Cos 0) 


PsCe) 



i£^l 

(i+7r J 


0 ( 2 . 66 ) 


(2.67) 

( 2 . 68 ) 


F (0) = Fg(0) + Fg(0) 


(2.69) 


^F (0)' stands for the contribution of gravity force and 

s 

'F„(0)' for the contribution of surface-tension force to the 
s ' 

condensate drain-off rate . 

If 'F(©)' of expression (2.69) now replaces' F(0)' of 
expression (2.18) in section (2.5.1), both formulations for 
the circular and the polar profile cases will be exactly alike. 
Therefore, the same solution method as adopted earlier for 
the circular profile is applied to solve the present problem. 

Some representative values of ' h' vs. obtained 

from the solution of the above problem are given in Table 2,3.' 
The effects of both gravity and stjrface -tens ion have been 
considered. The local heat transfer coefficient ' h^' as a 
function of ' 0' is obtained as : 

hg= [3.4363 - 2.2905© + 0.7086 0^-0.0779 x [k] 


(2.70) 
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where , 

'K' is given by expression (2.41) 

The average steam-side heat transfer coefficient 
is calculated to be : 

Sg = (0.7991) X [kJ (2.71) 

Percentage increase in ' in comparison with the 
corresponding circular profile is given by, 

'S " from expression (2.71 )4S' from expression(2.44) 

•/. increase = 2 ^xlOO 

'Eg from expression (2,44) 

= 9.42 y. (2.72) 

To study the isolated contribution of gravity towards 
augmentation in for the polar profile over 'E^' for the 

corresponding circular profile, the following procedure is 
adopted. 

In the expression (2.69) for 'F(O)', the surface- 
tension contribution 'F (©)' is set equal to zero. Then, adopting 
the same solution procedure, as is xased for the circular 
profile, the following results are obtained. Values of 'h^ vs. 

'Q' are obtained at various discrete locations along the 
periphery of the profile. Some representative values of "h" vs. 

'0' are presented in 'Table 2,4'. 

'h' is given by, 

s 



hs - [l .5879 - 0.5665 0.1925 ft^-0.0280 x [k] 

(2.73) 

' Eg' is evaluated to be : 

Eg = (0.7830) X [k] (2.74) 

/ increase in in comparison with the circular profile is 

calculated as : 

•/. increase = x 100 = 7,22 */ (2.75) 

‘/.increase in 'E ' due to the lone contribution of surface 
tension = (2.72)- (2.75) = 9.42 *X -7.22 ‘/. = 2.20 */. 

(2.76) 

From expression (2.72), it is observed that using the 
polar profile, ' r=a© ^ causes an increase in 'E/ by 9,42 */ 
over the corresponding circular profile. A heat transfer 
coefficient factor 'HTCF^ is now defined as : 

'S' for the polar profile from expression (2.71) 

HTCF= 

'E' for the circular profile from expression (2.44) 

(0.7991) X [k] 

= = 1 ,0942 (2.77) 

(0.728) X [k] 

From the definition of 'HTCF J it follows that this factor, 

when multiplied by 'S' for circul-nr profile, will give 'S' 

s s 

for the corresponding polar profile. 
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Therefore, for the profile 'r=a© 'can be 

represented as : 

E = (HTCF) X (0.728) x [kI (2.78) 

where, 'K' is given by expression (2.41) 

For the polar profile, under consideration, variation 

of 'h' with '©' is plotted in Fip,. (?..i). Values of these 
s 

variables of 'Table 2.3' and 'Table 2.4' are used, respectively, 
to graphically represent the contribution of 'gravity + surface- 
tension' and the contribution of 'gravity only ' to the heat 
transfer coefficient. Vari stion of' with 'w' for th.e 

corresponding circular profile, usin^;, ' 2.1 , la uiuo 

drawn in the same figure. A 'computer- programme No.o i-' given 
in 'Appendix-^,' which has been used to draw the above figure. 

A glance at Fig. (2.3) gives a feel of the relative contribution 
of 'gravity ' and 'surface-tension' to the heat transfer 
augmentation in the polar profile tubes in comparison with the 
corresponding circular tubes. 

Q ^ 

2 5 3 Profile, Represented by polar cu rve 'r=ae ', 

Reflected on both __s ides of..it3_yertlc^al^_ : 

Fig. (2. 4. a) shows the actual shape of the profile. 

The corresponding circular profile is also shown in Fig. (2.4.b). 
The method for drawing these figures Is the same as that for 
polar profile 'r=a0' (Section 2.5.2). since, the procedure of 
solution is exactly the same as discussed in Section (2.5.2), only 
the important results, associated with the use of this profile, 

are presented. 



13000.0 


4 



[Do z^/m] 


FIG. 2. 3 VARIATION OF ‘hj WITH '0' FOR THE POLAR PROFILE V =a0’ 
AND FOR THE CORRESPONDING CIRCULAR PROFILE. 
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AND THE CORRESPOND! 



dx = 

y2 ae^ dO 

(2.70) 

S = 

p 

150.6241 a 

( / 

RMF = 

0.0209 

(C-.B!) 

AF = 

0.8558 

(Z.8.’ j 

Ro = 

72 ae® 

(2.85) 

^ " 

iL e-2e 

2a^ 

• 

CD 

Fg(e) 

PjS 

= ^ (-Cos 0) 

^1 

(2.85) 

F3(e). 

°1 _ -20 

(2.86) 


2)Li^a 


'Table 2,5 ' gives some representative values of 'h' vs, 'Q \ 
considering both gravity and surface-tension. For this case, 
'hg is given by, 

hg = [j5.2555- 12.57'.5w + 5.7751 ^ j x [k] 

'S' is then calculated as ; 
s 

Eg = (1.U77) y. [■-.] ! '.• : 

'L increase in 'E' in <: .■n! irison wl Ui !.lu' uoru' - ■ itiU ; ■ 

s 

circular profile 

= 57.13 •/. (2.89) 

Heat transfer coefficient factor = HTCF = 1 .5713 


( 2 . 90 ) 
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O 


'Table 2,6' presents some representative values of 

hg vs , 0 j considering the effect of gravity only. For this 

case, 

[2.1628-0.5822 9+0.1510 ©^- 0.0197 ©^] x [k] 

(2.91) 

Sg= (0.9715) X [k] (2.92) 

*/, increase in ' over the corresponding circular profile 
= 23..03 */. (2.93) 

The surface-tension contribution towards the increase in 
'SJ = 57.13 i. - 33.03 y. = 24.10 */. (2.94) 

Fig. (2,5) shows the graphical representation of 'h' 

O 

vs. '©' for all the above cases. The same method, as adopted 
for polar profile, 'r=ae' in Section (2,5,2) is employed 
to draw this figure. 


2.5,4 Development of some generalised expressions : 

From the discussions in the Sections (2.5.1) to (2,5,3), 
the following generalised expressions are developed, which can 
be applied to any shape of the tube -prof ile . 


(a) The local heat transfer coefficient h : 


h = (HTCF) X (0.728) x 
s 


1 + 0.2 


Si^'^sat.-^w) 


h. 


fg 


(n-1) 


Pj g \ 


11/4 


X 


(2.95) 



36000 . 


OATIWAY 



AND FOR THE CORRESPONDING CIRCULAR PRORLE . 



in expression (2.95), 
are given as follov/irjg 


H'? 


dif f ereni 


oroiiii 


Profile 

Circular 

r=a© 


r=ae 


1 .0 

1 .0942 
1 .5713 



J 


( 2.96 


(t) ) The cross-sectional area ' k ' : 

cn 


A 

cn 


(AF) X 


''‘cn-^c 


(2.97 


In expression (2.97), (A ) ^ represents the 
cross-sectional area of the circular profile and 'AF' for 
different profiles are given as follov/ing ; 


Profile 


Circular 


' r=ae ' 


r=ae 


© / 


1 .0 

0.9758 

0.8558 


( 2 . 98 ) 


The above generalised relations (2.95 to 2.9S) will 
be used in the ensuing sections for the cost on t.t.i.iaa ti on 
analysis of the conaensing unit of a specific case. 
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Table 2.2 Values of at different locations *0^ along 
the periphery of the condenser-tube, 

( Circular profile ) 


SI .No . 

0 /~rtadian ~f 

^3 r.v/m2 

1 . 

1 .6022122 

7623.808 

2. 

1 .8849553 

7500.937 

5. 

2.1991143 

7386.562 

4. 

2.5132734 

7197.425 

5 . 

2.8274325 

6928.484 

b . 

3.141 5916 

6574.778 

7. 

3.4557507 

6125,256 

8. 

3.7699098 

5559 .1 54 

9. 

4.08^t0688 

4829.826 

10. 

4.3982276 

3797.853 

11 . 

4 .6809700 

1745.718; 



5S (a) 


Table 2.3 Values of at different locations ^9^ 
the periphery of the condenser-tube . 

( Profile, r=a9 . v/ith gravity 
and surface-tension ) ' 


bl ,Mo . 

0 /"Radian 7 

hg rvjm^ 

1 . 

1 .6022122 

13098.766 

2. 

1 .8849553 

8802.841 

3. 

2.1991143 

8234.439 

4. 

2.5132734 

7832.231 

3. 

2,8274325 

7407,571 

6. 

3 .141 591 0 

6')23.467 

7. 

3 .4557507 

6363 .61 4 

8. 

3.7699098 

5706 .693 

9. 

4 .0840&88- 

4908 . 249 

10. 

4.3982276 

3825.664 

11 . 

4,6809700 

1856.422 


along 


°C 7 
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Table 2. A Values of h at different locations 9 ^ along 

s 

the periphery of the condenser -tube , 



( j 

, ( p rof ile, r=a9 , 

with gravity only) 

SI .No . 

9 [~ Rad ian_7 


1 . 

1 .6022122 

rlwh .4 30 

2. 

1 .8849333 

H4<'i0.3R.9 

3. 

2.1991143 

! ii''.04 .8 i->5 

4. 

2.5132754 

7B54 . 378 

5. 

2.8274325 

lk5l 

6 . 

3.1A15916 

6949 .970 

7. 

3.4557507 

CD 

% 

vO 

8. 

3.7699098 

5719 .31 3 

9. 

4.0840688 

4911 .908 

10. 

4.3982276 

3824.526 

11 . 

4.6809700 

1746.612 
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Table 2.5 Values of at different locations ‘ 9 ’ along 
the periphery of the condenser-tube , 


SI .No . 

Q ’ 

( profile, r=.'ie , wi ih 

9 r~i(adian 7 

iV i ty 

iii'i surface- bens ion) 

hs rv./rn^ °C_7 

1 . 

1 .6022122 


3o3n4 ,306 

2. 

1 .8849553 


1 71 31 .704 

3. 

2.1991143 


12734.500 

4. 

2.5132754 


10636 .567 

5. 

2.8274325 


9334.939 

6 . 

3.141 5916 


8304.248 

7. 

3.4557507 


7335 .934 

8. 

3.7699098 


6342.675 

9. 

4.0840688 


5264.478 

10. 

4.3982276 


3974.910 

11 . 

4.6809700 


1860.889 
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Table 2,^ Values of ’ at different locations ^9 ’ along 
the periphery of the condenser-tube . 



- ^ w ^ 

( profile, r=ae , wi' 

th I'ravity otily) 

21 .ho . 

9 / hudian_7 

s 

1 . 

1 .6022122 

12364.641 

2. 

1 .8849553 

11674.965 

3. 

2.1991143 

10965.827 

4. 

2.5132734 

10178.106 

5. 

2.8274325 

9324.292 

A 

W # 

3 .1 4l 5916 

8414.611 

7. 

3.4557507 

7453.368 

8. 

3.7699098 

6433.805 

9. 

4 .0840688 

5323.377 

10. 

4.3982276 

3999.463 

11 . 

4.6809700 

1779.375 



CHm’I'H-III 


COST-BASED OPTIMISATION OF THE CONDENSER TUBE- 
DIMENSIONS FOR CIRCULAR AND POLAR PROFILE TUBES, 


3.1 INTRODUCTION : 


As has already been discussed, the film of condensate 
formed over the surface of the condenser tubes, accounts for 
the major resistance to heat flow in the process of film- 
condensation over horizontal tubes. The average unit 
conductance for heat transfer on the steam-side, in its 
most general form, associated with this process, is given 
by the expression (2.95), i.e. 


(HTCF) X (0.728) x 

O 


X 


1+0.2 — n,; (n-1) 


’'fg 

2 3 ^ 

pj g ^1 ^fg 




1/4 


(3.1) 


where , 

'HTCF ' is a heat .transfer coefficient factor, as 
explained in Chapter-II and given by expression (2.96) for 
different tube-prof iles . 


In expression (3.1), *dj stands for the outer 
diameter of the tubes if they are circular and for polar profile 
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-tubes, represents the outer diameter of a circular tube, 

which has been moulded to generate the polar profile. In 
either case, however, ' d^ stands for the outer diameter of a 
circular profile , 

Relation (3.1) indicates that the average steam-side 

heat transfer coefficient is inversely proportional to 

Reduction in 'd ' , therefore, results in an increase in 'S' . 

o s 

The water-side heat transfer coefficient is given 

by [ 4 ] , 

iV = ( ^) X (0.023) X (Re)°-® x (Pr)°-‘^ (3.2) 


'HMD' in expression (3,2) is the hydraulic mean 
diameter of the tube -profile and is given by the relation 
(2.58), i.e,. 


HMD 


4 X AF X 

cn c 

PER 


(3.3) 


where , 

'AF' is an area factor, which is defined in Chapter-II 
and is given by expression (2.98), for different tube -profiles . 

Since, ^^cn)c~ ^ ^i *= n d^^ the 

expression (3.3) for ^HMD ^ reduces to , 

HMD = (AF) X d^ 


(3.4) 
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"" v/here, 

'd/ represents the inner diameter of the circular- 
profile tubes and in case of the polar prof ile -tubes, ' d/ 
stands for the inside diameter of a corresponding circular 
profile, that has been moulded to generate the polar profile. 
In both cases, therefore, ' d^ represents the inner diameter 
of a circular profile. 


Using relations (3.4) and (3,2), the expression 


for 'h' becomes : 

W 


V ( 3?) ^ (0.023) X (Re)0-8 x (Pr)°-'* 

i 

(3.5) 


where 


Re = 


(p„ 1^^) X HMD 


u 


w 


or 


Re = (AF) X 


p V° d . 
^w w 1 


4 , 


w 


and 


Pr = 




(3.6) 


(3.7) 


The flow velocity of cooling water through the tubes 
in expression (3.6) is given by, 


V» = 

w A 


cn 


(3.8) 
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where , 

« 

^ is the volume flow rate of cooling water and 'A " 

or 

denotes the area of cross section of the tubes. ^ A ' is 

cn 

obtained from the relation (2.97), i.e., 

Kr^ = (AF) X (A ) 
cn ^ ^ cn^c 

or, 

= (AP) X ( 5 dj) (3.9) 


Using (3.9), the expression (5.8) for* reduces to , 




f n j 2 \ 


(3.10) 


Substituting (3.10) in (3.6), the expression for *Re* becomes 

(3.11) 


Re 




Substituting the expression (3,11) for *Re* and expression 
(3.7) for ^Pr^ in (3.5) and simplifying we get. 




k,, 

( — sp-g) X (0.023) X ( 


^ Pw ,0.A 

" ^w ^ ^ ^ ^ 


(3.12) 


Expression (3.12), shows that the water-side heat 
transfer coefficient, i.e., '^h^ is inversely proportional to 

i o 

(di) . For the same thickness of the tubes, a reduction in 
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the outer diameter requires a proportionate decrease in the 
inner diameter, which, in turn, results in an augmented 
water-side heat transfer coefficient. 

Out of the remaining thermal resistances in the 
path of heat transfer from steam to the coo I in/'; water, the 
major contribution is rendered by the foulin// of tubes , 
which is mainly a function of water quality and, hence, is 
independent of tube -dimens ions . 'Ihe tube-wall thermal resistance 
is very small in comparison with other resistances, encountered 
in the process. Thus, it can be safely neglected in the 
analysis, for all practical purposes. 

Hence, any reduction in the tube diameter results 
in an increased overall heat transfer coefficient. 

For a given condenser heat-duty, the flow rate of 

cooling water through the tubes of the condenser, before and 

after any change in diameter of the tubes, remains unaltered 

for a fixed temperature rise of cooling water . But in order to 

maintain the same flow rate with reduced inner diameter of tubes, 

requires a higher velocity of flow through the tubes. Since, 

the pressure drop in the tubes is proportional to the square of 

flow -velocity, the power required for pumping cooling water 

through the tubes is increased. The pressure drop through the 

tubes is given by [ , 

f L P 
2 X H I© 


AP = 


( 3 . 13 ) 
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Using expression (3.4) for 'HMD', the relation (3.13) 
for 'ap' becomes : 


AP 


f L p 
w '^w 

rx'(AP) X d. 


(3.14) 


Substituting relation (3.10) for ' in (3.14), the 

w 

expression for 'ap' reduces to , 


8 f L p 
(AF)^X n^x (d^)^ 


(3.15) 


It is evident from expression (3.15) that any 
decrease in the inner diameter of the condenser tubes results 
in a significantly high pressure drop through the tubes. This 
is because 'Ap' is inversely proportional to (d^) . Hence, 
the power requirement for pumping cooling water through the 
condenser tubes is greatly enhanced. 

To sum up, a decrease in the external diameter of the 
condenser tubes, resulting in an increase in the overall heat 
transfer coefficient, requires less surface area for causing 
a given heat transfer rate in the condenser. Hence, the amount 
of tube material, required for the pxxrpose is reduced, giving 
rise to a reduction in the first-cost of the condenser. But the 
consequent increase in the energy requirement for pumping entails 
an increase in the energy-cost vshich is a running cost, associated 
with the analysis. Against these two conflicting requirements, 
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namely, the decreased first-cost and tiie increased running 
cost, the optimum tube diameter based on minimum total cost 
has to be determined. Since, the energy cost is a running cost, 
it varies with time and a futuristic approach should be made to 
account for this changing cost element with time over the 
total service-life of the condenser, 'fhe 'present worth method 
is employed to perform this cost analysis in the ensuing 
sections of this chapter. 

The following three sections, (3.2) to Section (3.4) 
prepare the ground for carrying out the detailed cost analysis, 

3.2 DETERMINATION OF THE OVERALL^HEAT TRANSFER COE FFICIENT.'U.p 

AND THE THERMAL RESTS TANGE DUE TO CONDENSER TUBE- 

FOULING FOR CIRCULAR AND POLAR PROFILE TUBES : 

The following three factors contribute to the overall 
heat transfer coefficient, 

(a) Steam-side thermal resistance, represented by 'R^. 

(b) Water-side thermal resistance, denoted by 

and 

(c) Miscellaneous thermal resistance which is attributed 

to the fouling of the condenser tubes and is represented 

To determine 'R', 'h^' is first evaluated by 
using expression (3,1), i,e,. 
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S = (HTCF)x(0.728)x 


St ^Sat.-S) 

1+0.2 1: (n-1) 


h 


fg 


X 


V, 1 3 


1/4 


R' is then obtained with the help of the following 


relation 


s = 


H, 


(3.16) 


For obtaining ' R^ , ' h^' is first determined with 
the help of the relation (3.12), i.e,, 


1 k , 4p„ i 0.8 0.4 

^ t IF (0.023)x (-VT!— > ^ > 




/ 



is then evaluated, using the following expression : 




(3.17) 


where and 'A/ stand for the 'outer surface area' and 

'inner surface area' of liae condenser tubes, respectively. 
They are determined as following ; 


For a condenser tube having 
length = L m. 
outer diameter= d^ m. 

and inner diameter = d^ m. , 
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ivj is gjven by, 
arxd ' is given by, 

= K L (3.19) 


Using (3.18) and (3.19), the expression (3.17) reduces to. 




(3.20) 


The rate of heat transfer in the condenser is expressed as : 

Q = X p X C X (T^-T ) (3.21) 

w p^ 0 1 

where, 

is the volume flow rate of cooling water through 
the condenser-tubes . 

' t/ and ' tJ represent the inlet and outlet temperatures 
of cooling water for the condenser, respectively. 


^ Q ' and ' \]^ are related through the following 
expression : 

Q = (U^ A^) X LMTD (3.22) 

where , 

'LMTD' is the log mean temperature difference^ explained 


in Appendix-1 . 



(3.23) 


Using (3.22), we get, 

U = —p-S 

o X LMTD 

^he expression which relates the sum of all the 
thermal resistances and * U is given by [ 4] , 


U 


o 



where, 

E R = R + R + R 
s w m 


(3.24) 


(3.25) 


Using (3.24), we get the expression for ' ZR' as : 

ER = (3.26) 

o 

Since, ^ is already determined by using (3.23), 
ER^ can be evaluated from relation (3.26). ' R^^' is then 
obtained by rearranging the terms in expression (3.25) as : 


As all the terms of the right hand side of (3.27) are 
already evaluated, ' R^ can easily be determined. 

For circular profile tubes, using HTCF =1 ' and ^AF=1 ^ 
in expressions for ' and 'h^', respectively, the miscellaneous 
thermal resistance ' R^^ is found to be 0.7281 x °C/W 

when the total thermal resistance, i.e, ' ER' is obtained to 
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be 1.0187 X 10"^ ®C/W. Data for various parameters 

for the case study given in Appendix-1 have been used for 
the above calculations. 

3.3 DETERI4INATI0N OF A FUNCTIONAL RELATIONSHIP BETWEEN 
THfi tO&i OF production 0^ EN!er0Y AND TIME " 

'Appendix-10^ gives the cost of production of 
energy at ^Panki Thermal Power Station, Kanpur 'for the last 
seven consecutive years i.e. from the year 1 9 79 —80 to 
1985-86 . Considering the energy cost per unit as the 
dependent variable and year of observation as independent 
variable, a polynomial relation is obtained between the two 
variables. A ^computer programme No. 7', listed in 'Appendix -11 
has been employed to obtain this polynomial relationship by 
using the method of least squares . The polynomial, so 
obtained is given by, 

CEPU (t) = 31.7742 + 0.5632 t + 0,7124 t^ (3.28) 

where , 

'CEPU' represents the cost of production of energy 
per unit (KW . hr.) in Paise . 

't' stands for the time -element given by the year of 
observation column in Appendix-10. Substitution of 

7 in expression (3,28) yields the CEPU for the 

years 1979-80, 1980-81,, 1985-86, respectively. 
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The above functional relationaFiii.) is extrapolated to 
predict the energy cost per unit for the subsequent years. 


For example, substitution of t= 8,9,10, gives the 

^CEPU^ for years 1986-87, 1987-88, 1988-89 respectively 


This expression (9.28) will be employed in Section (3.5 
for the determination of the running cost-element, Involved 
in the cost analysis. 

3.4 THE PRESENT WORTH METHOD ; 

A brief description of the present worth method, as 
employed in the cost analysis is presented in this section. 

The total cost, involved in any cost analysis 
broadly comprises of the following two categories : 

(a) First cost, denoted by * FC . 

and 

t f 

(b) Running cost, represented by RC . 

Since 'RC' is a time -dependent cost element, its 
effect is entirely different than 'FC 'on the determination 
of the total cost. In any comparative cost analysis procedure, 
the present worth method acts as an elegant technique to 
account for the time -dependent and varying cost element, 
i.e., 'RC'. This method determines the present worth of the 
running cost, which is to be actually incurred at a future time 

The following example explains the discussion in 


detail : 
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Suppose "the "total investment on an equipment for its 
procurement is 'FC' , i.e., the fixed cost. 

Let the time -dependent running cost -element, 
associated with the operation of the equipment be denoted by 

t t it 

RC and 1 be the life of the equipment in years . 

^RC', being time-dependent, is denoted with a suffix. 
At the time of procurement of the equipment, i.e, at time, 

^ it it 

t = 0, RC is represented as RC^ , the magnitude of which 


will obviously be zero. 

At the end of the 1st, 2nd, 3rd,.. *^l*th years, 

RC is represented as RC.^ , ^RC^\ 


It is required to find the worth of the amounts 'RC.^', 


* * it 

RC 2 , RC^ , which are to be spent at the end of 

1,2, 1 years, respectively, I'rom the date of procure- 


ment of the equipment, at time t=0. In other words, the presen 
worth of the above amounts, which will be spent in a future 
period of time, is to be determined. 

Let the present worth of ^ RC^ be denoted as ' PW^' . 

This means that if " PV/.| amount of money is deposited in the 
bank at time t=0 , then at the end of one year, this amount 
' , alongwith its interest, will yield an amount ^RC.^. 

Initial amount = PW.j 

Suppose the annual rate of interest on bank -deposits 
is ^IR* in percentage. 
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Total interest, earned at the end of one year =PW^x 
Hence, total amount received from the bank at the end of 
one year = PW^ + PV/^ x = PW^ (1+ 

PW^ ( 1+ ^ ): RC^ 

RC. 

= .r-T^T 
( 1 + 

* TR ' ' * 

Replacing by X , the expression (3,29) reduces to. 


PW 


RC 


1 


1 ~ (i+ X) 


(3.30) 


Similarly, if ' PV/^ is "the present worth of ^^ 2 * 
at the end of two years, this amount alongwith the interest 
will yield an amount ^ 

Total amount that ^ PW^ will generate at the end 
of one year = PW^ ( 1+ (3.31) 

Total interest generated at the end of 

second year = PW 2 (1+ ^ ^ (3.32) 


Hence, total amount, which ' PW^ will yield at the end of 
second year = (3,31) + (3,32) 

= PW^ (1+ X)^ 


(3.33) 
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But this amount is the same as ' RC^' from the very 
djBfinition of present worth, 

RC^ = (1+X)^ 


or, 



RC, 


(1+X)‘ 


(3.34) 


Extending the same concept further, the following 
results can be obtained. 


PW3 


RC3 

(1+X)^ 


PW^ 


RC^ 

(1+X)^ 


A generalised expression for the determination of 
present worth can be obtained from the above discussion as 
following : 


If, 

RC^ = Running cost, incurred at the end of the 'i'th 
year . 

PW^ = The corresponding present worth. 


Then, is given as : 


RC. 

PW = K 

1 /4 . vNI 


d+x)- 


(3.35) 
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Hence, the present worth of the running cost, 
incurred over the entire life-span of the equipment is 
given by, 

1 

= 2 P¥ (3.36) 

i=1 

where , 

PW^ represents the present worth of the total running 

cost . 

The total cost ' TC\ then, is expressed as ; 

TC = FC + PWy (3.37) 

The total cost is calculated, using expression (3.37) 
for different equipments or different versions of the same 
equipment and that alternative, which leads to the minimum 
total cost is selected as the optimum choice. 

This, in brief, epitomises the present worth method 
for cost analysis . 

With the discussions in Sections (3,2) to (3,4) 
the stage is nov/ set to undertake the complete cost analysis . 
The details of this analysis are presented in the next section. 

3.5 DETAILS OF THE COST -OPTIMISATION ANALYSIS FOR CIRCULAR 
AND POLAR PROFILE TUBES Tn A CASS STuBy ; 

3.5,1 Circular Tubes ; 

The cost-optimisation analysis is carried out as 


following ; 
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In the case -study, the outer diameter of the condenser 
is 22mm. and thickness of the tubes is 1mm., i.e. 

d^ = outer diameter of tubes = 22 mm. 

t = thickness of tubes = 1 mm . 

The outer diameter of the tubes is then changed 
continuously and its effect on total cost of the condenser 
(fixed cost + running cost) is observed in order to obtain the 
optimum outer diameter for minimum total cost. The corresponding 
other optimum dimensions can then be determined . 

Hence, ^ d^' is a known quantity, initially. 

For a fixed thickness of tubes, the inside diameter 
* dj^^ is then given by , 

d^ = d^ -2t (3.38) 

The tube-wall surface temperature ’ T^^ has already 
been determined in Section (2,3) . With the knowledge of d^, 

'd/ and ' T,/ , the overall heat transfer coefficient'^ U/ is 
evaluated by following the same procedure as explained in 
Section (3.2). 

# • # 

The heat transfer rate Q in the condenser and the 
overall heat transfer coefficient ' Uj are related by the 
expression, 

Q = (U A ) X LMTD 
o o 


(3.39) 
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V.Taere, 

is the outer surface area of the condenser tubes 
and ' LMTD ' represents the log mean temperature difference. 

Q and *LMTD *are given in ' Appendix -1 1 


Hence, from expression (3.39), the outer surface area 
for heat transfer will be given by, 


^o "U X LM® 
o 


(3.40) 


But, as given in expression (3.18), 

A = 7c d L 
o o 

Hence, the length of the condenser tubes is given by. 


L = 


nuT 


(3.41) 


The volume of the tube -material is then given as : 


i (d^ - df ) X L (3.42) 

\Vhere, 

¥■ - Volume of tube material 

m 

The mass of tube-material used, is then expressed by. 


M = » X p„ (3.43) 

m m m 

Where, 

= Mass -density of the material for tubes 
M„ = Mass of tube -material 
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The total cost, incurred in obtaining the tube- 
material for the condenser, i .e . ' C^' is given by, 


C = 
m 


M X C 
m 


mtl 


(3.44) 


Vihere , 

*^mtl~ tubennaterial in Rs per kg. 


Since, the above analysis is carried out considering 
only two tubes (one for inlet and the other for the outlet 
of cooling water), the total cost of the tube-material is 
calculated by taking into account all the 13000 tubes of the 
condensing unit. Hence, the total material cost for the tubes 
is given by, 

TC^ = X (13000/2) (3.45) 

Where, 

TC = Total cost of the condenser tube ^nater ial . 
m 

The other element of cost, involved in the analysis 
is the energy cost for pumping cooling water through the 
condenser-tubes . The details of calculations, involved in 
determining this running cost element are as following : 

The pressure drop for water, circulating through the 
tubes, is determined by using expression (3.14), i.e., 

- ^ ^ ^ Pw 

2 x(AF)x d^ 
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in the above expression is taken as unity for 
circular profile tubes as given in expression (2.98) . 

The only unknown in the expression for *AP' is 'f'f 
which denotes the friction coefficient. If the flow of water 
through the tubes is laminar, no difficulty arises in 
calculating ' t \ because it is then given by the simple 
expression [ 6 ] , 

f = -I---- (5.^6) 

Actual calculations, however, show that the flow 

4 5 

Reynold’s number is very large (of the order of 10 to 10 ) 
and hence the flow is invariably turbulent. Several simplified 
correlations have been developed to obtain ' f for turbulent 
flow, but the most commonly used engineering practice is to 
use the coolebrook’s equation [6 ], which forms the basis for 
drawing the Moody’s diagram. This equation is given by, 

1- = - 2.0 log.. ( 4^ ^ (3.^7) 

^ 10 3.7 yj 

In the above expression 'e ' denotes the roughness 
coefficient for the tube -mate rial . Values for different 
materials have been given in Appendix -12 . 

Since, equation (3.47) is implicit in the variable 'f 

t § 

a and error ^ method is needed to solve it for f ♦ 
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The most logical way of solving this problem is by using 
the * Newton-Raphson's numerical technique % which is a very 
fast-converging method. 

In the present analysis, subroutine FRIXON of 
'computer programme No. 8', given in 'Appendix-13 ' calculates 
*f* by 'Newton-Raphson*s numerical technique'. 

Since all the parameters on the right hand side of 
the expression for 'ap' are now determined, the pressure drop 
'ap' through the tubes of the condenser can be calcxxlated. 

The power, required for pumping cooling water through 
the condenser tubes is then given as ; 

^w “ ^ 1 

lo 

Where, 

$ • t 

represents the volume flow rate of cooling water, 

if i 

given in Appendix-1 and 

* * stands for the overall efficiency of the pump, 

o 

which is "assumed on the basis of practically encountered 

( 

value to be 0.78 [ 4 ] . 

As has already been mentioned, the analysis is carried 
out for two tubes only. For the complete condenser, having 
13000 tubes, the total power requirement for pumping is 
given by, ! 

TP^ = (13000/2) [W ] 


( 3 . 49 ) 
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As the condenser is continuously running (24 hrs . a day 
the electrical energy, consumed for pumping per year is 
expressed by, 


TP X 24 X 365 

> _ 

'p ~ 1000 


KW. hr. 


(3.50) 


The cost of energy which is charged from a consumer 

towards the electricity bill is mainly a matter of Government poU 

and hence, it does not really reflect the actual cost of 

energy. The cost of production of energy, on the other hand, 

is more representative of its actual cost and therefore, is 

used in the present analysis. 'Appendix-10 ' gives the data for 

the cost of production of energy for the last seven consecutive 

c 

years from 1979-80 to 1985-86 , collected from the PANKI 
THERMAL POWER STATION, KANPUR . 

Suppose, the year for the procurement of a new 
condenser is considered as 1986-87. 

The approximate life of the unit is 25 years 
(Appendix 1 ) . 

Cost of energy, to be incurred at the end of 1986-87 
is given by, 

CEi = E X CEN^ 

* Jr 

Where, 

CEN^ = Cost of electrical energy per unit (KW.hr.) for 
the year 1586-87 (1st year) 
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and Ep has been determined from relation (3.50) 
dimilarly , for the subsequent years ^ 

CE^ = Ep X CEN^ 


CE 25 = Ep X CEN23 


In all these above expressions, CEN^ , CEN^ ^^^25 

are obtained by replacing ^t ^ in expression (3,28) by 

S»9, 32 respectively. The expression (3.28) is given as: 


CEPU (t) = 31 .7742 + 0.5632 t + 0.7124 

Using the ' present worth method ' , discussed in 
Section (3.4), the following calculations are made. 


With the help of relation (3.35), the present worth of 
the energy cost, which is to be incurred at the end of 1st year 
is given by, 


PWCE^ 


CE., 


where, 

PWCE = Present worth of energy cost 

= Rate of interest on bank -deposits, assumed 
as 0,1 . 

Suffix' 1 ' indicates the end of 1st year. 


X 
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Similarly, 

CEo 

PWCE^ = L_ 

^ {H-X)2 


PWCE25 = 


CE, 


25 


(1+X) 


therefore, the total cost; of energy, incurred over 
the entire life of the condenser, as represented by the 
pi'esent worth method is given by the following expression : 
25 

TCE = 2 PWCE /100 (3.51) 

^ i =1 ^ 

In the above expression, the factor *100' is used 
to convert the cost from 'paise' to 'Rs'. "^CEp' represents 
the total cost of energy for pvimping. 


The total cost, involved in the process is then obtained 
by adding the material cost and the energy cost and is 
represented by, 


TO = TC„ + TCE^ (3.52) 

m p 

where , 

TO » Total cost 

and ' TC^' is already calculated by using expression (3.45). 

A 'computer programme N 0 . 8 ', listed in 'Appendix-1 3 ' 
calculates the total cost for different tube diameters and the 



results are presented in 'Table 3.1' and 'Table 3.2'. In 
this tabular representation of the results of the cost analyst 
e=0.l5 mm. implies that the roughness coefficient of the 
condenser tube-material is assumed to be equal to that of 
galvanised iron and e= 0,046 mm. assumes this coefficient 
to be equal to that of ' commercial steel . 'These assumptions 
are made because of the absence of data for 's' for the 
condenser tube -material , which is ^Aluminium Brass. 

The optimum tube dimensions and related results of 
the cost analysis are then calculated by using the following 
procedure : 

At present, tubes of 22 mm. outer diameter are used 
in the condenser, taken for the case-study. Hence, all the 
parameters, given in Tables 3.1 and 3.2 are read against 
the row, denoted by 'do= 22\ These parameters are represented 
with suffix 'e' attached to them, which denotes existing . 

'when the total cost, denoted by ' CT' in the above 
tables, changes its nature from decreasing to increasing 
direction, this point of inflexion shows the optimum point. 
The optimum parameters of the cost analysts are then read 
against this point and suffix 'op' is attached to these 
parameters to denote optimum . 

Total aenrice life of the condenser =25 years approximately 

(Appendix-1 ) . ■ 
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for the various parameters is given as follov/ing : 
optimum outer diameter of condenser tubes 
Saving in electrical energy consumption per year 
Percentage Saving in electrical energy consumption 
Saving in energy cost over 25 years 
Annual saving in energy cost 
Percentage saving in energy cost 
Loss in material cost 
Percentage loss in material cost 
Saving in total cost over 25 years 
Saving in total cost per year 
Percentage saving in total cost 

The savings, indicated in the above parameters are for 
the optimum parameters over the existing parameters , 

'The expressions for evaluating the parameters (2) to 
(11) are given as : 


SEECY 

= (ENCYg-ENCY^p) x 6500* 

(3.53) 

PSEEC 

= (ENCY^-ENCY^„)/ENCY^ x 100 

6 op ^ 

(3.54) 

SECT 

= (CENTg- CENT^p) x 6500 

(3.55) 

ASEC 

= SECT/25 

(3.56) 


^ 'The factor 6500 is used due to the following reason, 

%e results of the cost analysis, presented in the "Tables 

5 1 to 3»6^ are on the basis of two tubes in the condenser. 
Hence, to take into account all 13000 tubes of the condenser, 
all parameters of the tables are to be multiplied by 
(13000/2) or 6500. 


Nomenclature 

(1 ) d = 

op 

(2) SEECY = 

(3) PSEEC = 

(4) SECT = 

(5) ASEC = 

(6) PSEC = 

(7) LMC 

(8) PLMC = 

(9) STCT = 

(10) STCY = 

(11) PSTC = 
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PSEC = 

(CENTg-CENT^p)/CENTg x 100 

(3.57) 

LMC 

(CMT^p-CMTg) X 6500 

(3.58) 

PLMC = 

(CMTop -CMT^)/CMTg x 100 

(3.59) 

STCT = 

(CTe -CT^p) X 6500 

(3.60) 

3TCY = 

STCT/25 

(3.61) 

PSTC = 

(Cl’e -CT^p)/CT^ X 100 

(3.62) 


In all the above expressions from (3,53) to (3.62), 
the magnitudes of the parameters of the right hand sides are 
re.-id from the' Tables 3.1 and 3.2.' These parameters are also 
explained in the nomenclature of these tables. The results, 
calculated by the above procedure are given as following : 

(a) circular profile, e = 0.15 mm [ Table 5.1 ] 

d^^ = 52 mm. 
op 

SEECY = 2887913 KW.hr., PSEEC = 99.^ ‘A , SECT = 263573000 Rs.A 
ASEC = 105429^2 tW- , PSEC = 99.4 */ , LMC = 4035070 Rs,/- 
PLMC = 56,43 */ , STCT=259488000 Rsy-,STCY= 10379556 Rs./- 

FoTC := 95.25 */, 

(b) circular profile, £== 0,046 mm [Table 3.2 ] 

= 51 mm. 
op 

SEECY = 2191.810 KV/.hr., PSEEC= 99.2 */. , SECT= 200041000 Rs./- 
ASEC = 8001669 Rs./- , PSEC== 99.2 */. , LMC = 3867565 Rs./- 
PLMC = 54.1 •/ , STCT=1 961 70000 Rsy-,STCY= 7846802 Rs./- 

PSTC = 93.9 '/. . 
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Polar Profile Tubes : 

In the expressions (3.1), (3.12) and (3.14), which 
ave been used in the cost analysis for circular profile tubes, 

' and ' HTCF' are replaced by the following values for the 
olar profiles; 

Polar profile AF HTCF 

r = aO 0.9758 1^0942 

r = ae® 0.8558 1 .5713 

The above values of ^AF ' and 'hTCF' have been already 
valuated in Section (2.5,4). 

The rest of the procedure is entirely the same as 
as been adopted for the circular profile tubes in Section (3.5) . 
ithout repeating the procedure, therefore, only the final 
esults are presented, 

(a) Profile, V=a0 ^ , £=>0.15 mm. [ Table 3.3] 

d^_ = 52 Hun, 
op 

SEECY» 2991925 KW.hr., PSEEC= 99.4 */. , SECT » 273066000 Rs./- 

ASEC = 10922657 Rs./-, PSEC » 99.4 •/. , LMC « 3907670 Rs./- 

FLMC = 54.7 */. , STCT = 269054000Rsy-, STCY=1 07621 90 Rs./- 

PSTC « 95.5 •/. 

(b) Profile, ^ r=*a0 \ 0.046 mm, [ Table 3.4] 

d = 51 mm. 

^op 

SEECY “ 2271101 KW.hr , PSEEC='99.2 */ » SECT = 207278000 Rs./ 

ASEC » 8291137 Rs./- , PSEC =99.2 •/. , LMC = 3740880 Rs./- 
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PLMC = 52.3 •/. , STCT = 203433000 Rs./-, STCY = 8137342 Rs./- 

PSTC = 94.2 */ 

[ Table 3.5 ] 

SECT = 339860000 Rs./- 
LMC = 4260295 Rs./- 
Rs./-, STCY=1 342401 5 Rs^- 

[ Table 3.6 ] 

SEECY = 2829085 KW.hr., PSEEC = 99.3 '/. , SECT = 258204000 Rs./- 

ASEC = 10328177 Rs./-, PSEC = 99.3 */. , LMC =4090905 Rs./- 

PLMC = 60.8 i. , STCT = 254113000 Rs./-, STCY =10164541 Rs./- 
PSTC = 95.26 */. 

Another aspect of the results of cost analysis 
[Tables 3.1 to 3.6 ], which is used to find the saving in 
tube -Biate rial by using polar profile tubes in the condenser 
over the corresponding circular profile tubes, is presented 
as following : 

'The condenser taken for analysis in the case -study, uses 
circular tubes of 22 nm. outer diameter. If these tubes are 


(c) Profile, ' r=ae® e= o.l5 mm. 

% = 54 mm. 

SEECY = 3723774 KW.hr., PSEC = 99,5 */. 
ASEC = 13594427 Rs./-, PSEC = 99,5 •/ , 
PLMC = 63.4 y. , STCT =335600000 

PSTC = 96.3 y. 

(d) Profile, r=ae® e = 0.046 mm. 
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changed to the polar profiles, then there is a reduction 
in the tube -material, required for effecting the same heat 
transfer rate in the condenser. This causes a decrease in 
the space -requirement of the condenser and also, there is a 
reduction in the first cost. 

The actual magnitudes of the above results are 
obtained for polar profiles ' r= a© ' and ' r=ae® * respectively. 

PROFILE ^ r = a© ^ 

Saving in material = ( MTM from Table (3.1) - MTM 

from Table (3.3) x 6500 

= (10-9.855)x 6500 = 9^2.5 kg. 

where, ' MTM' represents the mass of tube-material 
and is read against the row, given by ' D0=22' 
in the tables. 

*/ saving in material = */• 

Saving in first cost of the condenser= Saving in material x cost 

of material 

= 942.5 X 110 

= 103675 Rs./- 

The savings, indicated above are for the polar profile^ r— a© 
over the corresponding circular profile. 
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profile V = ae^ ' 

Saving in material = ( mTM from Table (3.1)- MTM from 

Table (3.5) x 6500 
= (10-9.AO6) X 6500 = 3861 kg. 

Where, MTM denotes the mass of tube -material and is read 
against the row, given by ' DO = 22' in the tables. 

Saving in material = x 100 = 5 ,9^ '/. 

Saving in first cost of the condenser = Saving in material 

X cost of material 
= 3861 X 110 

= 424710 Rs./- 

'Ihe savings, indicated above are for the polar profile, 

'r = ae®'over the corresponding circular profile. 
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Nomenclature for Tables (3.1 U n 


DO 

FF 

j^NCY 

CENT 

GMT 

CT 


Outer diameter of tubes £~rm._7 

Friction factor j/~dimensionless_7 

Energy, consumed per year 7~Kw.hr._7 

Total energy cost over the entire 
service life of the condenser /~^s ._7 

Total cost of tube material /~f^s._7 

Total cost, involved, for the entire 
service life of the condenser /” Us .7 


MTM : Mass of tube material K'.s, ._7 

L : Length of tubes M 

e ; Roughness coefficient of tube-material [ nran.] 


Note : This analysis is made on the basis of one set of two 

tubes (one each for entrance and exit of cooling water) 
of the condenser. There are 6500 such sets in one 
complete condensing unit. Results for the complete 
condensing unit, therefore, are obtained by multiplying 
these results by 6500. 



Table 3,2 Results of cost analysis 
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CHAPTER-IV 


DISCUSSION OF RESULTS 


A numerical technique to solve the problem of 
determining the condensation heat transfer coefficient for 
any shape of tube -profile in a condenser has been developed 
in this thesis. The heat transfer augmentation effects, 
obtained by using the polar profile tubes, given by polar cxirves 
'r=a0' and 'r^ae® in comparison viih circxiLar profile tubes 
in a condenser, are presented. The coirre spending reduction in 
the first cost of the equipment is also evaluated. The optimum 


condenser -tube dimensions are determined on the basis of 
minimum total cost of the equipment. The condensing unit of the 
110 MW system of the ' Panki Thermal Power Station, Kanpur has 
been taken as a case-study for the above analysis . Solutions for 
this specific case have been presented in the form of tables and 
figures (Tables 2.1 to 2.5, Tables 3.1 to 3.6, figures 2.2 to 
2.5). Important results on the basis of these tabularly 
represented solutions have been obtained in Chapters II and III. 
In this Chapter, a detailed discussion of these results is 

presented as following : 


For ocwenlence of representation, tne polar profile of 
cor^ienser tuOes. given Oy 'r-aO' la represented as ' PPCT-1 ' 
and for the cu^ye ' r=ae« L » is denoted hy ' PKT-a 
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The use of 'PPCT-1 ' in place of circular tubes in 
the condenser results in an increment of 9.42*/. in the heat 
transfer rate /"Chapter -I I_7. 7,22 */. of this increase is due 
to the contribution of gravity force and the rest 2,20 /. is 
because of the surface-tension force. Hence, the effect of 
gravity is predominant in causing the heat transfer augmentation. 
The use of PPCT-2 causes an increase of 57,13 */. in heat 
transfer, out of which, 33.03 */. contribution is due to gravity 
force and 24.10 */. is the effect of surface-tension force. 
Therefore/ both effects have almost even contributions for 
causing the heat transfer augmentation in case of PPCT-2 * . 

As has been pointed out in Section (3,2), the thermal- 

resistance across the heat-flow path in the condenser due to 

tube -fouling is found to be Q.7281 x 1(5^ m^ ^C/W when the 

—3 

total thermal -resistance is 1 .0187 x 10 . Hence, fouling of 

condenser tubes, accounts for 71,4 */. of the total thermal 
resistance, which is quite large. This warrants regular cleaning 
and maintenance of the tubes to keep. this value of thermal — 
resistance equal to the usually accepted value in commercial 
practice ( of the order of 0.176 x 10 ^ m^ °C/¥ jkj ). This 
will ensure a high heat transfer in the condenser and hence, 
a smaller size condenser is required for a given heat^^uty. 

From the results, presented in Section (3.5.1), it is 
obsenred that using' pPTC-1 ' and ' Pi>TC-2 respectively, instead 
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of the circular tubes in the condenser, result in a reduction 
in condenser-space requirement of 1 ,45 V. and 3.94 */. respectively. 
This, in addition to the discussion of the previous paragraph, 
establishes the fact that whever, in an industrial system 
using a condenser, space is a primary concern, * PPTC -2 * may be 
used in place of circular tubes in the condenser. 

The most important results, however, are those, obtained 
from the cost-optimisation analysis in Section (3.5), Consider- 
able saving in total cost of the condenser has been seen to 
accrue from the use of optimum tube -dimens ions over the existing 
dimensions, ^or example, changing the outer diameter of the 
condenser tubes from 22mm. to the optimum diameter, which is 
found to be 52 mm. for e=o.l5 mm /“Section 3.5__7, will cause 
an annual reduction in total cost, amounting to Rs .10379556/- 
f or the entire service life of the condenser , 

Moreover, using the optimum diameter of the condenser 
tubes, mentioned above (52 mm), also results in an annual 
saving in electrical energy, amounting to 2887913 KW.hr. over 
the existing tube-diameter (22mm.). This is a fabulous saving 
in electrical energy, which is a very precious commodity, 
specifically in the wake of the present energy-crisis . Hence, 
adoption of optimum tube dimensions in the condensers of all the 
thermal power plants of our country will result in an 
incredibly large saving in energy, which can be diverted to 
other energy -starved areas. 


many 
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Similar results of the cost-optimisation 
analysis for 'pPTC-1 ' and 'PPTC-2' have been presented in 
Section (3.5.1). 

An apparent demerit in the use of proposed polar 
profile tubes in place of circular tubes in the condenser 
seems to be that the manufacturing cost of these unconventional - 
shaped tubes may be very high and this may cause a large 
increase in the first cost of the equipment. But in all 
probability, if the manufacturing processes for these tubes are 
standardized, the cost, involved in their manufacturing will 
be comparable with that of the corresponding circular tubes 
and therefore, not much increase in the first cost of the 
condenser is expected. 
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APPENDIX -1 


OPiAPv TlNG DATA OF THE CONu ENAKH, 
UNIT AT THE PANKI THlPGMyVl. PnVJ!-:H 


U/iED F.-'U 
STaI'ION, 


AHA 110MW . 
KANPUK - 


(a) Condenser pressure (Kgf/cm^) : 

Max . 

Min . 

Normal 

(Vacuum) : 

0.95 

0.80 

0 .85 

(.Absolute) : 

0.103 

0.233 

0 .183 


(b) iAoolint', water temperature 

at the inlet to the condenser : 32°C = 

(c) Cooling water temperature 

at the outlet from the condenser: 4o°C = 

(d) Tube size : 

(i) Outer diameter = 22mm. = 

(ii) Inner diameter = 20 mm.= d^ 

(iii) Length = 7.5m. = 

(e) Material of Construction 

of the condenser tubes : Aluminium Brass . 


(f) Volume flow rate of cooling 
water through the condenser 

(g) No. of hours, the condenser 
operates in a day 

(h) Approximate life of the 


15AOO m^/hr. 


24 hrs.' 


condenser. 


25 yrs 



( ij I'ic). ol Lubes in the condenser 
(j) Arrangement oC tubes 

) JiVerUi'^e number of tubes per 
vertical row. 


: 13,000 

: 2-tube-pass, shell 
and tube arrangement . 

: 40 


U) Vanlght of one tube of 
7.3 m. length. 

(in) dost of pipe material 


5 Kg. 

Rs . 110 per Kg . 


OA’l'A, DERIVED FROM THOSE, SPECIFIED ABOVE :- 


(n) Saturation temperature of 
steam in the condenser at 

normal pressure. : T , . = 57.76*^0 

S 4 A c 


(o) Latent heat of evaporation of 
steam at the normal pressure 


h 


2364064 /j/Kg. 


(p) Length of tube, to be considered 

for evaluating pressure drop : ]., = 15 in. 

(q) The volume flow rate of cooling water through 


one set of tubes is calculated as follows : 


Total volume flow rate of cooling water 


3 



V/;it:e.r' liov^s in through one set oC tubes and flov/s out 
oi. the condenser through another set of tubes. Hence, 
only half of the total number of tubes are to be 
considered for calculating the flow rate through 
one-pass . 

Volume flow rate of cooling water through one-pass 
(consisting of one tube from each set, above) 

_ _ 15AOO 

= 2.3692307 m^/hr. 

2.3692307 
- “ 36^50 

= 6.58119 xIO”"^ m^/sec. 


(r) fhe mass -density of the material of condenser tubes 
is calculated as follows : 

( ) 

A tube of 7.5m. length with dQ=22mm. 
and = 20mm.’ weighs 5 Kg. 

Volume of the above tube 

= 5 (d^- d|)x10“^ X 7.5 
-4 3 

= 4.948 X 10 m 

Mass -density of the tube material is then given by, 


ra 


5 Zp Kg./m^ 

4.948 x10' 

10105 Kg./ra 
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(y) I'he LMTD (log mean temperature difference) is 
evaluated as follows : 


LMfl) 


(T 


sat 




log. ( , 


T - T 
sat. o ^ 


__p 

^sat. -^i 


(57.76-40) - (57«76-32) 

,57.76-40 ^ 

l°Se ^7. 76-31 ^ 


- 8 

_ /' 17 . 76 ^ 

loge (25 .IE’ 

= 21 .51 °C 


(t) 


Q = 


Heat Transfer rate through I5m. 
condenser = ^ 


length of the 




-T^)=1 000x4 


X 4.184 ( 40-32) 
» 21890 ^ 
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appendix -2 

'.rHEHMOPHY S ICAL PROPERTIES OF CONDENSATE (WATER) /~ 3_7 


Nomenclature 


T 


P 


C 


P 


k 

4 


Temperature /~°C _7 
Mass -density /“Kg./m ^_7 
Specific heat /~J/Kg. °C _7 

Thermal conductivity /"W/m °C _7 

Dynamic Viscosity _/“N.s/m ^_7 


Interpretation ; 

To use the last column of the data (shown on 
the next page in a tabular form), the following procedure 
is adopted. 

Let the value, read from this column for h 
be 792.4 . 

Then, the actual value of h is 

792.4 X 10“^ N.s/m^ 


. . . .Contd. . . 
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Ap]ietKi:i,x-2 (,contd ,) 


■r 

P 

c 

k 

(I X 10 

30 .0 

995.70 

41 76 .00 

0.615 

792.4 

32.0 

995.10 

4175.60 

0.619 

763.4 

33 .0 

994.10 

41 75 .00 

0.624 

719.8 

37 .0 

993.34 

41 75 .00 

0.628 

695.1 

ho .0 

992.20 

41 75 .00 

0.633 

658.0 

42.0 

991 .40 

41 75.40 

0.636 

636 .9 

45 .0 

990.20 

41 75 .00 

0 .640 

605.1 

47.0 

989.36 

41 76.80 

0.643 

585.1 

50.0 

988.10 

4178.00 

0.647 

555.1 

52.0 

987.04 

41 78.96 

0.649 

540.8 

55 .0 

985 .46 

4180.40 

0.652 

519.4 

37.0 

984.40 

4181 .36 

0.654 

505.1 

60 .0 

982.82 

41 82 .80 

0.657 

483.7 


E 
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APPENDIX -3 




if\w 1 

r 

(.•('Mf 'llTFfV pKOGRAMMf: No. 4 

001 ’' 1 1 

t 

■ • - 

00000 

i- 


00100 

i: 

PRObKAMMh to tVALUATP THK. FUNCTIONAL REL AT I ONHH 1 P 

00000 

c 

OF THERMAL CONDUCT I VITY (K) , DYNAMIC VISCOSII Y (MU) , 

00 500 

i 

Sf^ECIFIC HE At (CP) AND MASS DENSITY (RHO) WITH THE. 

00400 

c 

VARIATION OF TEMPERATURE USING THE 'METHOD OF LEAS f SQUARES'. 

004 1 0 

(' 

♦#*#**#** ^ Sr ■¥ ¥;**■* 

00000 


DIMENSION S( 10 , 10 ) ,X< 100 ) ,B(5) ,A(5) ,Y( 100 ) ,CY< 100 ) , PD (100) 

000 10 

c 

SSSS-SSSrSrSr-¥fWSrSr^*-SS:S-SS;S-SS**S*-¥i-****-¥r** SSS-SSSSStSS-SSSS 

00520 

c 

NOMENCLATURE: 

00^00 

t,; 

M=^-NUHB£R OK DATA POINTS 

00700 

c 

N*ORDER OF POLYNOMIAL TO BE FITTED THROUGH THE DAI A POINTS 

0 OB 00 

c 

X^aNDEPENDENT VARIABLE (TEMP. ) 

00 ^?O 0 

c 

Y«T‘)EPENDENT VAR I ABLE (RHO , MU , K OR CP) 

01000 

i, 

CYwDEPENDENT VARIABLE CALCULATED FROM THE POL YNQMI AL -F 1 ! . 

0 1 \ 00 

r; 

PD*PERCEN 1 AGE DEV I AT I ON 

0 1 200 


SUM^CUMULATl V 6 'PD' . 

0 1 1 i*) 

( 

#4f ###.«• ##4t 4*- 4*### ^f****##*-)<-***##4f *##***« *******■« s SrS-S^^¥--ih*¥rSS-SSS-SrS 

01220 


WRITE (49, 51 ) 

01250 

*1 

T ORMAT (///////////////) 

0 1 ;^. 0 O 


m 1 3 

0 1 400 


m- 2 

01407’ 

t 

###+« 4< 4f 4f .«■#* .)^ ******* .K-* *** 4«- 4<- **** **. .X- .♦f. ** +• 

0 1 *>00 

0 

THF FOLLOWING FILE GIVES INPUT DATA TO THE PROGRAMME 

01507' 

f 

+« >!■****** 41 ********+( ***^< >r**** ****** *^» *^« ****** >^* ****■> **** 

O 1 7*00 


RFAD{24,*) (X (1) ,Y(I) ,1^1 ,M) 

OlctO'* 

i 

H X M * M .V * H'**** •»< M**** ***■»* ^* ^ * X* **** *#***^1 **** K*^***********^^^* 

01 *00 

i; 

LVAiAtATiari OF THE ELEMENTS OF THE MATRIX ’S' AND MATRIX 

01 ,\V* 

» 

-^<*4^2l^*X*X*■*f•KX'|lt*^<******************^<*********^<*XX***^**********^‘*» 

0 1 MOO 


Ot jMY W0 . 

i* 1 ’*'*^***2 


tiOMXN :0. 

07000 


01 IMXSY.i'0. 

O. UH' 

t 

HI fMXLY -0. 

02200 


HUMX-0. 

07 200 


SUMXi;h''0. 

02 400 


r,HMXCB-'i:0. 

02l *V10 


MOMXFR"»0. 

0:v,v’io 

i 

SUMXFV«0. 

0 2 nu* 

I 

rutMXsx'v.0. 

0 2MOO 


DO 10 I* KM 

0. **00 


14 IMV- SlIMV LT ( 1 ) 

o *.000 


FU IM X V • Hi )M X V 4- X ( n -M Y ( I > 

0 ' 1 00 


HHMXLV vBUMXViV^' ( X ( I ) **2) ^Y ( I ) 

0 '200 

( 

SI IMXLV "250MXI:Y4- ( X ( I ) *Y ( I ) 

0 


SMMX- 14IMX'* X ( n 

0 '4 00 


F.;t JMXSR 'SUMXSR < X ( 1) 

0 ’MOi 


tiltMXLH’‘'BllMXLB'»-X ( 1 ) **3 

0 ■-,-.00 


Si IMXF R'".*SUMXFR"» X ( I) ^<-*4 

0 • '02. 

1 

HUMXt V SIIMXI V-+X(I>*^5 

0 ,,i02 

1 

‘SUM X S X "Bl IM X B X ■+• X ( I) * ^6 

0 "'00 

10 

LON f INUF 

04 020 


F» ( t , 1) '*M 

04 1 00 


5(1 ,2)^-BUMX 

04 200 


Sn ,5)‘^BUMXSR 

04 ".OO 

i, 

Bit ,4)‘.^-9UMXLB 

04 400 


S ( / , 1 ) «SUMX 

04 *.00 


iU2/2) wSUMXBR 

04 • 00 


S ( 2 % 5) •JsBUMXCH 

0.1 /0O 

i 

5(2,4) «SUMXFR 

04!40O 


S ( 5» 1 > IMXSR 

047 00 


S(5/2)"*SUH5(0H 

014000 


S (5, 5) ^SUMXFR 
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0510(3 

C 

S (3, 4 ) -BUnXhv 



05200 

C 

S (4, 1 ) -SUMXCB 



05300 

C 

S<4,2)=SUMXFR 



05400 

c 

S(4,3)=SUMXFV 



05500 

c 

S{4,4)=SUMXSX 



05600 


B ( 1 ) =SUMY 



05700 


B<2) =SUMXY 



05800 


B(3)=SUMXSY 



05900 

c 

B (4) =surixc:Y 



06000 


N=N+1 



06100 


WRITE (21 ,*) < <S ( 1 , J ) , J”1 ,M> ,1 = 1 

,N) 

, (B ( I ) , l«t ,N> 

06200 


CLOSE (UNI T=21 , DEV I CE= ' DBf ’ ' > 



06300 


READ (21 < (S( I ,J> , J= 1 , T = i , 

N) , 

<B ( I ) , 1=1 ,N) 

06309 

c 

^******#***^<.*»***** *******■«■** 



06400 

c; 

INVERSION OF THE MATRIX ‘S'. 



06409 

c 

********•)<•******-»************* 



06500 


1=1 



06600 


NX-N+l 



06700 


NY=2'»^N 



06800 


DO 80 wl=NX,NY 



06900 


S < I , J ) =1 . 



07000 





07100 

so 

CONTINUE 



07200 


L=1 



07 300 


K = 2 



07400 

110 

XM=S(L,L) 



07500 


DO 140 J=L,NY 



07600 


S(L,J>=S(l ,J) /XM 



07700 

140 

CONI INLIE 



07800 


DO 190 I=K,N 



07900 


X1=S(I ,L) 



08000 


DO 191 J=L,NY 



08100 


9(1,J)=S(1,3)”'5(L,J)-*X1 



08200 

191 

CONTINUE 



08300 

190 

CONTINUE 



08400 





08500 


T=K+1 



08600 


IF (L~N) 110,1 10,230 



0B700 

230 

L-~N 



08800 

2 55 

LZ=L“1 



08900 


DO 290 K= 1 , 1. 2 



09000 


I=L~K 



09100 


Y 1 =S (1,1.) 



09200 


DO 291 J=L,NY 



09300 


S(I ,J)=S(1 ,0)*Y1 



09400 

291 

CONTINUE 



09500 

290 

CONTINUE 



09600 


l.=L-i 



097 00 
09900 

320 

] r ( L- 1 ) 320 , 320 , 235 

WRITE (23, *) C (S ( I ,3) , J=NX ,NY) , 

1 = 1 

,N) 

09900 

1 0000 


C LOSE (UNI 1 "-23,01-: VI CE= 'DSK' ' > 
READ ( 23 , * ) ( ( S ( I , 3 > , 0 =NX , NY ) , I 

= 1, 

N> 
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10009 C 
10100 C 
10109 C 
10200 
10?.00 
10400 
10500 

10600 71 

10700 70 

10 70<=.' C 
10800 C 
10809 C 
1 0900 

1 1000 
11100 
1 1 1 09 L 
1 1 200 C 
1 1 :.>00 ( 

1 1 509 c: 

1 1 A 00 

1 1500 
1 1600 
1 1 700 
1 1 800 

1 1 900 28 

1 2000 
1 2 1 00 
1 2200 

12’S00 2/ 

12400 
1 2500 


MULTIPLICATION OF THE INVERTED MATRIX OF 'S' WITH MATRIX 'B'. 

***'#*^f**^**************************#****-K-'M-**-*******-»(-*-if-**4>e-***-)f 

DO 70 1=1 ,N 
A(1)=0. 

DO 71 a=NX,NY 
A*< I ) =A{ 1) -»-S ( 1 
CONTINUE 
CUNT INUE 

4«.***4t* )t^4^.*^(*^****4(.4(.*4(.**^*******^t *4( ***^*4t**** ■# * 

THE ROW MATRIX 'A' GIVES THE COEFF . S OF THE FITTED POLYNOMIAL. 

WRI TE (48, *) (A < I) , 1 = 1 , (NY-N) ) 

CLOSE (UNI r=48,DEVICE='DSK ' ) 

READ (48,*) (A( 1 ) , l^l , (NY-N> ) 

*4( >f.***4t 4^ **•>«•*■»»• #■# ** 

CALCULATION OF THE DEVIATION BETWEEN TFIF ACTUAL VALUE AND 
CAi. Cl HATED W\1 UL Of THE DEPENDENT VAR 1 API t*. 

BUh=0. 

no 27 D=1,M 

CY(v3)=0. 

DCJ 28 I = 1 , N 

CV ( J ) =CY (J) + (A(I)*<X(v))**(l-‘n>) 

CONTINUE 

PD ( J ) ^ ( ( CY ( J ) -Y ( J > ) / Y ( D M * 1 00. 

SUM=SUM^PD < J ) 

Wf a T E ( 49 , * ) X ( J ) , Y ( J ) , CY ( vT ) , PD ( .1 ) , SI »M 

CONTINUE 

STOP 

END 


. t ype 


\ 
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APPENDIX -h 


POLYNOMIAL 

FirmFnTS' 


RELA TIONS BETWEEN THE THERMOPHYSICAL 
mi COm'NsATE ( water] 'A ll) T'MPE'FUTIMe 


1. P = 1000.922 - (0.1757813 x 10“”^ )T -(0 .5875752x10“^)T^ 

2. Cp= 4196.031 - (1 .010742)1’ +(0.1684570 x 10"'^)T^ 

3. k = 0 .5383759+ (0 .3154755 x 10"^)T-(0 .191 3682 x 10”^)T^ 

4. ui = 1371 .895 - (23.90112)T +(0.1526718)T^ 


Note : Values of thermophysical properties, listed 
in''Appendix-2' and ^computer programme No.l' 
given in Append ix- 3 ^ have been used to develop 
the above corelations. 



Ill 


APPENDIX -5 

COMPUTER PROGRAMME No, 2 


00100 C 
00200 C 
00300 C 
00400 C 
00500 C 
00600 C 
00700 C 
00800 C 
00900 C 
01000 C 
01100 C 
0 1 200 C 

0 1 300 C 
01400 C 
01500 C 
01600 C 
01700 C 
01800 C 
01900 C 
02000 C 
02100 C 
02200 C 
02300 C 
02400 C 
02500 C 
02600 C 
02700 C 
02800 
02900 
03000 
03100 
03200 
03210 
03215 U 
03300 
03400 
03500 
03600 
03700 
03B00 
03V00 
04000 
04100 C 
04200 C 
04300 C 
04400 C 
04500 C 
04600 

04700 28 

04800 

04900 

05000 

05100 

05200 

05300 

05400 

05500 


PROGRAMME FOR OBTAINING THE TUBE-WALL SURFACE TEMPERATURE 
BY ITERATIONS USING ’IF' LOOP. 

♦ ♦ ****#***##*#**•>»■*♦##♦♦*♦ *♦*•»*** *********** ■‘^ * ^ ■'*'**■■*'■**■* “^ * "** 
IL«INITIAL LENGTH 

K=THERMAL CONDUCTIVITY OF LIQUID FILM 
MU^DYNAMIC VISCOSITY OF LIQUID FILM 
CP»SPECIFIC HEAT OF LIQUID FILM 
RHO=MASS DENSITY OF LIQUID FILM 
HFG«LATENT HEAT OF VAPOURISATION 

HFGSC«'HFG' CONSIDERING SUBCOOLING OF THE LIQUID FILM 
DQ«OUTeR DIAMETER OF TUBE 
AO«OUTER SURFACE AREA OF THE TUBE 

HO«HEAT TRANSFER COEFF. OF CONDENSATION r^r^KvrscrMeATP ■ 

CCFSC=CHEN'S CORRECTION FACTOR FOR SUBCOOLING OF CONDENSATE 
^IvjX»AVERAGE no. of TUBES IN A VERTICAL ROW 
TSAT«SATURATION TEMP. 

TW*WALL TEMP. OF THE TUBE 

TREF«REFERENCE TEMP. r^oKintTMocp 

QOR«0RIGINAL HEAT TRANSFER RATE IN THE CONDENSER 

QNEW«NEW HEAT TRANSFER RATE AFTER ITERATION 
ADG«ACCELERATI0N due to GRAVITY 
PCD«PERCENTAGE DEVIATION BETWEEN QOR AND 

FRHO FK,FMU,FCP ARE FUNCTIONS DEPICTING THE VARIATION 
W^TH TEMP OF ' RHO * , ' K ' , 'MU AND 'CP' RESPECTIVELY 

FRHO (xUl 000. 922-0. 1757813E-01*X-0. 3875732E-02* ( X**2) 

FCP < X ) «4 1 96 . 03 1 1 . 0 1 0742 *X +0 . 1 6e4570E-0 1 * ( X ) 

FK <X) »0"538i759-»-0. 3154755E-02.X-0. nr^L) 

FMU(X)«<1371.e95-23.90112*X^0.15267l8*(X**2))/(10.**6) 

WRITE(40,U) 

FORMAT (//////////////''^ 

HFG*2364064.0 

TSAT=57.76 

PI»3. 14159 

ADG®*9« 81 

DO»*0.022 

IL«15.0 

ao^pi^do-^il 

QOR»21090.0 w^.**^***********************'*^*'^**^* 

************************** given and the -if- loop 

^TAPTS^^”E^eRHfNATXcPITEPION IS THAT 'PCD- IS LESS THAN 
OR EQUAL TO 0.B1 '’^^^f^*****- 

TW*54. 0 

XPiEP*TW'+-0. 25* (TSAT-TW) 

MU“FMU (TREF) 

K«FK<TR£F) 

CP»FCP (TREF> 

HreSC-^HFsf/c (3. /a. > * (CP) * <TSAT-TW) > 

ant- 40. O ,.rp,(TSAT-TW>*(ANT-1.0)>/HFS 

HOn=”RH0**2) * (HFGSC) * (ADO) * (K**3) 
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05600 

05700 

©5800 

05900 

06000 

06100 

06200 

06300 

06400 29 

06500 9 

06600 
06700 
06800 
06900 

.k/i 


H0F2= (MU) * (TSAT-TW) * (ANT) * (DO) 

HOF« ( HOF 1 /H0F2 > . 25 

<0. 728) *(CCFSC)*(HOF) 

QNEW=: (HO) * (AO) # (TSAT~TW) 

PCD= (ABS (QQR-QNEW) /QOR) *100. 0 
IF (PCD.LE. .01 )60 TO 29 
TW=TW-f 0.001 
GO TO 28 

WRITE (40,9) TW,HO,QOR,QNEW,K,MU,RHO,CP,CCFSC,HFB8C 

FORMAT (15X, ' TW« ' FB. 5/ / iSX , * H0« 'F8.3//15X, ' QQR- 'F8.2//13X, 'QNEW* 

lF8,2//i5X, 'K«'F8.6//15X, ' MU« ‘ E13. 7/ / 15X , ‘RHO» 'FB.3//15X, ' CP- ' 

1 F8 . 3/ / 1 5X , ' CCFSC« ' F8 . 6/ / 1 5X , ‘ HFGSC- ' F8. 0 > 

STOP 

END 
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0t300S' t; 

0m 1 c 
00011 (' 

00 1 00 C 

00?00 ( ' 

00300 r 

00*100 r. 

00^00 c 

00600 L 

00700 C 

00H00 C 

00 ^ 0 V^ (■ 

01000 I 

0 1 1 00 I 

0 1 '.: 00 i 

01 *00 0 

0 1 1 00 1 

0 1 m0 I ' 

0 1 <^.00 I 

01700 ( 

01000 I 

0 1 '*’ 01 - ( 

0/000 t,. 

0. 100 ( 

0/000 

0? '.00 I 

00*100 i 

00! ^10' C 

00600 i; 

00 ? 00 * [ 

00900 C 

0OV0O ( 

03000 i: 

0':’100 I 

0 • .'*00 C 

00 '.00 L. 

0 0400 (,: 

0‘*O:.00 c; 

00.600 1 

0 0 700 f. 

0 *.000 1 ’ 

0 "O'^OO t 

04000 ( 

04 too 
04 000 
04 ;.00 
04 400 
04* .00 
04/.00 
04 '00- 
04 000 


COMPIIVER PROGRAMME No.O- 


<*** if *#■»(■ K ^♦-^<.************ 


PKDBRAMMfc ID OBTAIN THE HEA'I TRANGf-lrR CfJtFMCIEN] 
AT VARIOUS DISCRETE LOCATIONS ALONG THE PERIPHERY 
Of- THT- CONDENSER TUBES OF ClRCDLAfv PROf ILt . 




NOMENL'L.A 1 UUi : - 


RMF=RADTIIS MUI 11 plying factor. 

f'iAD-ROiL' n I! . (If Ulf f JRCUI 04.; f'fvUf h | . 

A--CONE}T Ar4 1 IN 1 Hf F’Ul Nf* 1 x Pf . I ' >1 i 1 HU Ul U ii f'li.p 

f ( I HE Up) . 

I REF ~:KE f- F. fvl', NCE T f: Mf '{ R A 1 Ufvl- . 

RHCll '-lit- NSl r . 01 I I UUI D. 

T L - T HE RMAl i .ONI ) I ( ' H V I 1 t' Ul I 1 (.4 1 1 u . 

MUI -im^AH 1 1 VISIUSIIV Of I lot I ID. 

S I Gl SI IRf Ai ’I T F N‘-; I UN OF I Hf I 11 M M 1 u 

HfU-lAUNi HTmI Of VmPUUR 1 U M < AJ IntitiwUiiiU ft.f* .* K-l . 

HF GSC- ’ HF t-: l ’l K IS I L)F R I NG 1 HI F f f f ( I Of ' .0| M ( lUl 1 r-F . Of (HI 


(:t)NDLfT..P. I f . 


Cl f Bt -.LHl 0 ‘ t.(il4.;n llliM f r,( lUf.' f UR' SUDM K 4 lUf.. 

AN 1 -64Vf-:RF,(iI llllHDFfv Of JUDIP PH'* '.4-'R(!(^4 R( U.J Ol- 


IHF C-ONDfNF.! 1 >. 


r SA r-BAl I IRO 1 I Of 1 I FMF’ERA 1 1 IRF . 

1 Ul- T FMPEfV'.UIfrF Of IHf WAl I bllf4A(f Ul I Hf lHf4^ . 
ADH-- ACCEl-ERFi I 1 UN DUf I O UPAV 1 1 . 

i;UNH«-(,.( »NS I AN I 111 Hf APPI HI) fU U II f t INi , 11 UN H 
WHICH CUNIRI HUIFb 111 I HP GRAVITY F.HHl UN UK 
lUNDhNFlAlt DF-AIN-UIH f.‘Atf'. 

CONSOl -CUtvlSI r.Nl lU HF AF'Pl IhD lU IHf- FUNt 1 I UN FD 
WHICH CUN IK I HU IPS 1 U THE SURf ACE - 1 F.NB U IM EO EL (. I 
OH THF: CONI )E N'oA T E DRA ] N-OU I RA TF . 

X-ANGl.E JN RADIANS, MF.ASORED FROM THE UPPFf- 
(tJENFlRA I R 1 Uf I HE POLAR PROF 1 LE . 


NU-No. OF UDSt.RVAl IONS. 

I N«VAU )l:, 01 1 HE I N TEGRAT 1 ON . 

(; p SK/FU.ISED IN EVALUATING THF INTEGRA!,. 

H r (:=:HF.A'( T RANSF ER CGEFF I C 1 E.NT . 

DEI 1 A'-HOUNDaRY layer THICKNESS OF THE LONDFNbAlE 


FORMED OVER I HE TUBE-SURF ACE . 

NOIE:*- ALL VARIAHL.ES ARE MF:ABI,IRH1) IN 'Si' ONllS. 




real F , 10, in, I I ,MUI 
D I ME NS I ON Y 100) , H T C (100) 

F AH. 

F*B ( X 1 ='0. (CONE' H (0. 5) * (F-'Xr* < *'/ ) > 

I ( < V > s- o:i,iNi'>) f CL IB ( X ) ) 

FT. C X ) -T 17 ( > ) I H t 

FI c ) ... <rn 1 - ^ 

^ ,, XH*n (Y 1 #FAC<> 
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ir,.t 

062»e' 

0".. ';0D 
05400 
0 t '.500 
05 A 00 

0 r '.:'’00 

0 tie 00 
0 f:. 9»0 
0 vbCU 10 
0 ot 00 

0 V -.200 
.*00 
0*',‘U*>0 
0 ^* - 00 - 
06. *>00 
tV'- ‘00 

0 a :'«00 

iv ‘>•00 
0 7 000 1. 

0 M 00 » 

0 ’ ;‘00 I 

0 •* .00 f 

0 '400 
0 '.00 
0 . . ,S 00 

0 V '00 
0 /B00 

0 ' 700 , 

0>;U^0l'^ 

0H100 

09000 

00r,00 

0*0400 2 H 

gi0S00 

0't ^iS00 

0f'!..'00 

0s'n-n't0 

0 1 .•<,.‘00 
0*7000 

0 ‘ 7 t 00 10 

0''7.‘0t1 
ili'/ ’100 

09400 
0*7500 
iVt'.600 
09 700 
0« ‘1*100 
0 'r. 4 l 0 
1 0000 
» 0 i 00 
1 0200 
U 5 ' t ' 1 


WK 1 I H ( , :'. > 

FORMAT < / ."' //////// / 

K’MI- = 1.0 
RAD=.01 1 
A-- (KAD) * (RMK ) 

RH0l.=9B8. 1638 

M =0. 654 r;.? 

MLIL = ( 5 - 1 80925 1 / ( 1 0 . * )^4 ) 

EUGL---.07 4 
HF0=.:‘'364064 . 0 
HI- (:i8L-=236V5V.l . 0 
CCFSC=1 . 048099 
A^4^ =40. 0 
I 8A'f =57. 76 
I W=54 . 2B 
l> ] rSA1 - 1 W 
Al)Cv='7. 81 

Cl INC 1 = ( RHl H. HI- GSC ( ADO ) 

r'l INOV •- t 3 .)<*■( M ) -I* ( Mill ) ^* nrn l A ) 

Cf}NB=C:ON(it '( (1H8.? 
l :Or48 I 0: I 81 ( a i ♦> ( Kt ll U ) +• I MM 4i« 1 
CI.INS f 2-- ^ < M ) * ( A# ^ ^^0 *» l MMI ' « ' D M 

CON81 «C;nM''M 1 M l INS! ? 

EVAIIIAIUUI 01 IHl INrt'Olv'A' 'IN H V MO. INi. f.niF’StifJ 
Rl )1 .f: F OR: NMMC- 1^ 1 1 f\\ I N MM IRA I 1 1 IN . 

X 1 M f. ; -M? V 
. 602/1 2.' 

NO=0 

M™> ( X2"' X n , . . 

1-2 

S MM' o: n 4 !• M’ 2 
S2«0. 

B4=K O* 1 +H) 

I N “ I ?> 1 4 . +t 8 A 1 * ( H /' . ) 

I F ( AH8 { ( I N” 1 0 ) / I N ) . I F. . . 00 M Gl 1 10 29 

R.2~'S2 + 84 
S4=0. 

X'-XMMI/V. 

DO 10 4^1 , 1 

8ArS44R (X ) 

Xi’X4-H 
(,„c:ir4 1 1 r4UF 
H'“*H / 2 . 

1 -2*)- .1 
ID* IN 

I N'"’ C M 1 4 ..M fe 4 4 . S 4 ) * ( H / . 1 

I,n iA:--( I UN)*^* <(?».::■&) ir;') n .(Tim-.*.'.'?) 

HTt (n = 'Kl 'I.)L-.I -tA)*(CCFSf;5/(ANr**0.7S> 

V ( n^'-Hic< I ' 

WP n £ ( 24 , ■H 1 V ( u 
N0=ND+1 

WP ri R ( 7’ , M ND , X 2 , H T C (I ) 

K c IPMA2I ( / 20X , 1 2 25/ . • 2 , f.-X 2' 'i' - ^ 



115 


1CJ300 C 
10400 C 
10500 C 
10600 C 
10700 C 

10800 c: 

1 0900 C 
11000 
11100 
11200 

11300 59 

1 1400 
11500 C 
11600 C 
11700 C 
1 1 800 C 
11900 C 
12000 
12100 
12200 
1 2300 C 
P 

12400 i: 
12500 C 
1 2600 C, 
12700 l: 

1 2900 C 
1 2900 t 
1 3000 C 
13100 C 
1 >200 C 
1*5 500 I 

13400 c: 

1 :5I-i00 C 
25600 C 
i 3700 C 
1 3800 C 
13900 r... 

14000 C 
14100 C 
14200 C 
14300 C 
14400 C 
1 4500 9 

14510 C 




THE following FILE STORES THE SOLUTION, GIVING THE 
HEAT TRANSFER COEFFICIENTS AT DIFFERENT DISCRETE 
LOCA*T IONS, SPECIFIED BY 'X2' AND THE TOTAL No. OF 
SUCH REPIT IT I VE CALCULATIONS ARE PERFORMED BY THE 
THE 'IF* LOOP THAT FOLLOWS THE 'WRITE* STATEMENT. 




IF (X2.GT.4.680970)GO TO 59 
X2«X2-»'0. 0314159 
BO TO 39 


CLOSE < UN I T-24 , DEV I CE- ' DSK.. ' ) 
READ(24,->t) (Y(I) ,1 = 1,99) 




THE FOLLOWING FILE STORES THE HEAT TRANSFER COEFFICIENl 
DAT A, TO BE USED IN THE PROGRAMME FOR FITTING A POLYNOMI- 
AL THROUGH THESE DATA POINTS. 


■»»-****************************#*#*****#**#**-j(-*****^«-*****}j-**** 


WRITE (24,*) (Yd) ,1 = 1 ,99) 

STOP 

END 


**************^4t*********4(.********-»^********-«**->t-»t****** ********** 


NESSARY CHANGES, WHICH ARE TO BE MADE IN THE ABOVE PROGRAMME 
TO USE THE SAME FOR SOME OTHER POLAR PROFILES, ARE AS 
FOLLOWS* — 

(«<); FOR THE PROFILE ' R=A^< THETA ' : 


1. THE EXPRESSION FOR 'RMF* 

IN 

LINE 

No. 

'4900 ' 

IS 

CHANGED TO 'RMF*0. 3021 ' , 

2. THE EXPRESSION FOR *FA(X)' 

IN 

LINE 

No. 

• 4 300 ' 

IS 

CHANGED TO ' F A ( X ) =SQRT ( 1 . 0+X**2 ) 
3. THE EXPRESSION FOR 'FBCX)' 

IN 

LINE 

No. 

' 4400 * 

IS 


CHANGED TO * F B ( X ) = (CONST )*(< X**3+4 . *X )/(< 1 .* X**2 ) **3 > > . 

(b): FOR THE PROFILE ,R*A*EXP (THETA) ' i 

1. THE EXPRESSION FOR 'RMF ' IN LINE No. '4900* IS 
CHANGED TO * RMF=0 . 02085 * . 

2. THE EXPRESSION FOR ■FA(X)* IN LINE Nd. ‘4300’ IS 
CHANGED TO ' FA ( X ) =SQRT (2. > *EXP ( X ) ' . 

3. THE EXPRESSION FOR ‘FBCX)* IN LINE No. '4400 IS 
CHANGED TO ‘FB (X) = (CONST) * (0.5.) * (EXP (-2. *X) > ' . 

TO STUDY THE ISOLATED EFFECT OF GRAVITY ON 
THE HEAT TRANSFER COEFFICIENT AUGMENT AT ION , IN BOTH 
THE ABOVE CASES '(a)* AND 'b*,THE EXPRESSION FDR 
FB(X) IS MULTIPLIED BY 'O.0- TO NULLIFY THE CON- 
TRIBUTION OF SURFACE-TENSION. 

*******#4^.M-********************************^************'*^***’**^ 


.type y 
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c; 

Qm 1 0 c 

00011 c: 

00049 c 

00050 t. 
00052 C 
00054 C 

00056 C 

00057 t 
00100 

00 1 0 1 f 

00102 t: 
00102 C. 

00104 r: 

00 1 05 L 

00106 C 

0010/ c; 

00 1 0B C 
00 1 09 C* 

00110 C 

00111 C 

00112 r 
00200 
00300 
00400 
00300 
00600 2 
00700 

00 7 08 C; 
00710 C 
007 1 2 C 
00/14 C 
00800 
00900 

01000 
01 100 
01 200 

0 i 2 1 0 c: 

01300 
01400 
01500 

0 1 600 

01 ;00 
0 1 800 
01810 t: 
01820 C 
01 V00 
02000 
02100 
02200 
0*2300 
02:310 I. 
024 00 


COMPUTER PRQGRAHMP No. 4 


*#*■>«•******•«•■)*•* •« •H>t'**'******-IH<-*-»f-****'lH**->H*!»»*-»H» 

RRD6RAMME TD DETERMINE IHE FUNCTIONAL REI.AIIQN 
BETWEEN THE STEAM-SIDE HEAT TRANSFER COEFFICIE" 

Nt HS’ AND ANSLfc ' T HETA ' , ME ASURED FROM T Hf: UP- 
PER GENERATRIX, USING THE 'MEVHOD OF LEAST SUUARES . 

*■)<*•>< #***«■<♦•* 

DIMENSION S< 10, 10) , X (100) ,B (5) , A (5) » Y < 100) , CY < 100 ) ,PD< 100) 

4i. 4(. * ^ -K 4< -X ** *<<•■♦*■* # -J*- *K *■»<**#*■«*•«•* ^ *****•**•* -X ^ > 

NOMENCLATURt : - - 

M«No. OF DATA POINIS.- 
X« 'THETA' ,MeASUR£D IN RADIAN. 

Y« 'HS ACTUALLY OBTAINED. 

N«QRDER OF THE POLYNOMIAL TO BE PITTED. 

CY«'HS' , CALCULATED FROM THE POLYNOMIAL 
RELAriON. 

PCD-PERCENfAGE DEVIATION BETWEEN 'Y' AND 
'CY ' . 

BUM=^CUMUL AT I UE ' PCD ' . 


M«99 

X(l) = 1.602'2122 
1)0 2 0 = 1, M 

X<J>=X(1)+(J~1)*^0-0314159) 

CONTINUE 
N=:3 

#*^«' ****•*•«+ 

FOLLOWING FILE READS THE VALUES OF 
DETERMINED FROM ANOTHER PROGRAMME. 

44**4,.*4t*******-»^********- 

READ<23,*) (Yd) 

SUMY»0. 


n 4< ■»***■«•****• 


SUMXY»*0. 

SUMXSY=0. 

SUMXQY=0. 

SUMXFYa©. 


BUMX»0- 

SUMXSRW0- 

SUMXCB«0- 

SUMXFR»0. 

SUMXF V«0. 

SUMXSX*0- 


SUMXSN»0. 
SUMXET=0- 
DO 10 l«liN 


SUMV»SUMY*+-Y (1) 

SUMXV^SUMXY+X (D^Vd ) 

SUMXS Y=SUMXSV-»- ^ ^ ^ ^ J 

SUMXFY»SUMXFV+ (X d ) ) *Y d ) 
SUMX=*=SUMX + Xd ) 



021'J00 

02600 

02700 

02800 

02900 

02910 

02920 

03000 

03009 

03010 
03012 
0301^ 
030 i 5 
03100 
03200 
03300 
03400 
03410 
03500 
03600 
03700 
03800 
03810 
03900 
04000 
04100 
04200 
04210 
04 300 
04400 
04500 
04600 
04610 
04620 
04630 
04640 
04650 
04660 
04700 
04B00 
04900 
05000 
050 1 0 
05 1 00 
05 1 09 
05110 
051 1 2 
051 14 
05'200 
05300 
05400 
0!.»409 
0! *4 1 0 
0t .41 1 
01.500 
05.600 
05700 
05800 
01..900 
06000 
06 1 00 
06200 
0630&"' 
06400 


5 

t: 

10 

L 

C 

c 

c 

c 


c 


I 




5 

c 

c 

c 

t; 


r: 


i. 

c 

c 


5 

r 

V 


£91 

1 10 


SUhXSR«SUMXSR+X < 1 > 
SUhXCB=SUMXC8+X il )** 5 
SUHXFR=SUMXFK^-X ( I ) 
SIIMXFV=SUMXFV4X ( 1) 4(*5 
SUhXSX=SUMXSX-+-X 1 1 ^ *#6 
SUMXSN=SUMXSN4X ( I ) 
SUMXET=SLJMXer + X ( I ) **8 
CONTINUE 


#*•1^ ***-* ^t #***<♦** It* -It I. *#■*.* * 

•S' REPRESENTS THE COEFFICIENT MAT RI X ,GENF K‘A U-l) FROM 
THE SIMULTANEOUS LINEAR’ tDllAnONS,Oin AKsll I) Kih: LTUVi- 
NG THE LEAST SQUARE ME 1 HOD. 

K ************ #**!** It -M ***«***#****-*f#***#****K»i«***'*^«* 


S ( 1 , 1 > =:M 

S(1 ,2)=SUMX 

S(1,3)«SUMXSR 

S(i,4>--=SUMXr.B 

S( 1 ,5) =SUMXFR 

S<2, 1 ) =-SUMX 

S(2,2)=SUMXSR 

S(2,3>=^^BUMXCB 

S(2,4)«SUMXFR 

S(225)=SUMXFV 

S(3, 1 ) =SUMXSR 

?~i(3,2>s=SUMXCE< 

S(3,3)«SUMXFR 

S<3,4)=SUMXFV 

S(“5,S>«SUMXSX 

B(4,1)«SUMXCB 

S(4,2)»SUMXFR 

B(4,3>«*SUMXFV 

S<4,4)«SUMXSX 

S<4,5)=SUMXSN 

S<5, 1>*SUMXFR 

S(5,2)«SUMXFV 

S(5,3)=SUMXSX 

S<5,4>*BUMXSN 

S<5,S)«SUMXE1 

EK 1) =:SUMY 

B <2) =SUMXY 

£i<(:5)*=SUMX6Y 

B<4)«SUMXCY 

B(S)=SUMXFY 


********** * 4 

* H^FOLLOriN/ K a-F. SIOREB THE ELEMENTS Of- THE 
CtlEFFlClENT MAIRIX 'S' 

5a ,N> , J=1 W'l. 1 (FUI ) , 1 = 1 .N) 

CLOSE (UNn=21.t)eVICE='DSK-) 

E:EHLI (21 , *> 1 IH ( I , J ) ,0- 1 ,N) , I - 1 ,N) , 

^ ^ M ************* **^ * ‘>^'»*"^-»^*-^'-*< •“ 


, 1=^1 ,N) 


} riVE FVS 1 ON OF MAT R I X ‘ S ' = 


I -t 

NV»N+ 1 
NY-L-^'N 

1.0 80 J = NX,NV 
5 0 = 

5-0 4 1 
f.UNT ir4UE 
( =0 

xri=s<i. iL) 
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06500 


DO 140 J=L,NY 

06600 


S(L ,J)=S<L,a) /XM 

06700 

140 

CUNT INLIE 

06800 


DG 190 I=K,N 

06900 


X 1 =S(I , 1 .) 

07000 


DO 191 3 =L,NY 

07100 


S(l,J)=S<i ,vT)-B(l 

07200 

191 

CONT I NUE 

07300 

190 

CONTINUE 

07400 


L^L^-l 

07500 


F.=--K.+ l 

07600 


IF <L-N> 110 , 1 10 , 2 ' 

07700 

230 

L=N 

07800 

235 

L 2 =L -1 

07900 


DO 290 - 1,1 7 

08000 


It=c-K 

08100 


Y 1 =S <1 , 1 ) 

08200 


DG 291 J=L,NV 

08300 


S < J , v 3 ) =S ( 1 , J ) -S F 

08400 

29 J 

CONTINUE 

0 B!t 00 

290 

CONT INUE 

08600 


L=l-l 


08 7 W 
09 B 00 
0 B 900 
09000 

09008 

09009 
0901 1 

09 1 00 
09200 
09 300 
09400 
09500 
09600 
09 700 
09800 

09808 

09809 
09010 
0981 1 
09900 
0990 1 
0990 V 

0990 

09904 
0990 !) 
10000 
10100 
10200 
10300 
10400 
1 0!-.00 
1 0600 
10/00 
1 0800 
10900 
1 1000 
11100 


'•>20 


71 

70 


t : 

L 

C 

{ 

c 

c 

c 

t 


I » XI 


B . 
^ * ■» 


IF <L -1 ) 320 , 31 ‘ 0 , 23 t. 

WR I TE < 43 , < S (I , J ) , 4 «NX , N r ) , I - 1 , H > 

Cl.DSE (UNn=43,nEr.VlCF=-- 'DBl ' > 

READ( 43 ,^) < (S(l ,.V-‘NX,NV) , I* 1 ,N) 

MUL run ICAUON of inverted matrix of S with MAfKIX 

DO 70 1 = 1 , N 
AU )= 0 . 

DO 71 J«NX,NY 

A<I)=A(n^’S<l 

CONTINUE 

CONTINUE 

wRi TE <44 <A ( 1 ) , 1 - 1 1 ^ 

hAlRIX -A' STORES THE VALUES OF THE COEFF- I C I tN I S OF 

2Lr°*i*°*i**F**”****™°^^^^^ 

SUM=0. 

txi o=i,M 

CY <T>= 0 . 

rv ( j*) =CY un + (A < I ) • ( X ( j) ** ( i-n ) ' 

CUNT INUb 

s rop 

Efin 


. t V P € 
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COMPUTER PROGRAMME No*5 


00100 0 
00200 C 
00300 C 
00400 C 
00500 C 
00600 C 
00700 C 
00800 C 
00900 C 
01000 C 
01100 C 
01200 C 
01300 C 
01400 
01500 C 
01600 
01700 
01800 
01900 
02000 
02100 
02200 
02300 C 
02400 
02500 C 
02600 
02700 
02800 
02900 C 
03000 9 

03100 C 
03200 C 
03300 C 
03400 C 
03S00 
03600 
03700 
03800 C 
03900 C 
04000 C 
04100 C 
04200 
04300 
04400 
04500 
04600 C 
04700 99 

04800 C 
04900 C 
05000 C 
05100 C 
05200 
05300 
05400 
05S00 C 
0S600 C 
05700 C 


* **♦***■«•*****#***** -tt *•*<•*•♦(•** *#■!<■♦ **#***♦■>»•#***♦*•»♦■ ■«■*##**#*#■♦<■ * 

PROGRAMME TO PLOT THE POLAR PROFILE OF THE CONDENSER 
TUBES BY USING 'INTERACTIVE GRAPHICS LIBRARY < IGL ) ' 
SUBROUTINES (FOR PROFILE ' R*A*THETA ' ) . 

4*- ************* *****■»<■***•»•>♦•♦****** -N- #■**♦***♦*•»•**♦* ^^^ ★ '><'** ■*■ 


MAIN PROGRAMME TO DRAW THE. P.I.R.C. PROFILE. 
NOMENCLATURE; 

RAD^RADIUS OF THE CIRCULAR PROFILE, THAT HAS BEEN 
MOULDED INTO THE POLAR PROFILE OF THE FORM 
■ R™A#THETA ' OK ' R=A*EXP ( THETA) ' . 

THETA«ANGLE, MEASURED FROM THE UPPER GENARATRIX. 

P, I.R.C.*PR0GRESSIVELY INCREASING RADIUS OF CURVATURE. 


DIMENSION R(25) ,THETA(25) 

INITIALISATION OF THE GRAPHICS TERMINAL. 

CALL GRSTRT (4010,5) 

CALL MILLIM 

CALL VWPORT(0. ,100. ,0. ,100.) 

CALL WINDOWCO. ,100. ,0. , 100. ) 

TH£TA(1>«90. 

RAD»11.0 
A=0. 3021*RAD 
As::0.02085-«'RAD 
R(1)«A*(3. 14159/2. ) 

R < 1 > «A*EXP <3. 14159/2. > 

DO 9 J«2,21 

’I HETA ( J ) =THETA (D-fCJ-D^dSO. /20. ) 
r(J)o(A)*( <3. 14159/2- ) + (J~l)*(3. 14159/20. ) ) 

R < J ) » <A) *£XP ( (3. 14159/2. ) + ( J*“l > * (3. 14159/20. ) ) 

CONT X NUE ***********4f******»^***************^^* 

TH^FOLLOWINB FILE STORES THE SETS OF 
TO DRAW THE P.I.ft.C. PROFILE(LEFT HALF SIDE). 

I°22*,;*L*u****..*****************‘**************** 

WRITEt23,*) (R<I ) ,THETA(I> ,1-1,21) 

CL0SE<UNIT=23,DEVICE='DSK-) 

READ<20,*) <R<I) ,THETA(I) ,1 = 1 ,-^l> *#**#**.»*«***#* 

p^niRF TO DRAW THE P.I.R.C. PROFILECLEFT HALF SIDE). 
CALL FIQURE<R, THETA) 

?hEtLjI=360.*180.-1THEIA(1)h(J-1)MI80./2B.)) 

3 14lU/2. ) + (0-l)*'T. 14157/20. ) ) 

CONT I NUE « »i* 

™rro"*o:r»rFrurii;;;;^ 

THE P. ^ ’*^^^^^^^^*i.*5***»** »•.****•***'**■**'■**'*** ******* 

WR1TE(30,*) <R(J> ,THETA(J) ,J=1 

CLOS£(UNIT=30,DEVICE» DS) ) 

read( 30.*) <F!<J> ’■''•^^''^‘'iLdC.***************’' •*•******* 

********************* jpr,purU4(.- 'FlGURf DRAWS THf 
the FOLLOWING CALL TO 

P.I.R.C. PROFILECRIBHT HALF SIDE). 
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05800 ^ 

05900 
06000 ^ 
06100 
06200 C 
06300 
06400 
06500 
06600 
06700 ^ 

06800 C 
06900 C 
07000 
07100 
07200 
07300 
07400 
07500 
07600 
07700 

07800 '5»0 

07900 
08000 
08005 C 
08010 C 

08110 c: 

08210 C 
08310 
08410 C 
0B510 
086 1 0 C 

08615 C 

, t y \y t \ '' C 


*■«•*■«•>«•**•«• •«■**•»■«•■«•*■«■***** ***4t*4^*^l***** *#*#*#**#♦**« •I* *-«•■»* 

CALL FIGURE <R, THETA) 

* * *** ^ #****#* -ifr ******** *** * 

DRAWING THE CIRCULAR PROFILE. 

***** **•<*■* 4*' •♦I- ***■«■*#**■♦( * 

CALL MOVE (80.0,80.0) 

CALL ARC (RAD,0. 0,360.0) 

CALL GRSTQP 
END 

***************************************************** 

THE FOLLOWING SUBROUTINE 'FIGURE* IS USED IN THE 
****************************************************■<» 
SUBROUTINE F I 6 URE ( R , THETA) 

DIMENSION R(25) ,THETA(25) 

CALL PIVOT(50. 0,85.0) 

CALL MPOLARCRCl) ,THETA(1) ) 

DO 90 J=2,21 
CALL SMOOTH 
CALL GRAIN (0.0) 

CALL DPOLAR(R<J) ,THETA(J) ) 

CONTINUE 

RETURN 

END 

********** ************************************** 

THE FOLLOWING CHANGES ARE TO BE MADE TO DRAW THE 
POLAR PROFILE ' R=A*EXP (THETA) ' BY USING THE SAME 
PROGRAMME. 

DETACH THE LABEL C*, THAT IS, THE COMMENT 
STATEMENT FROM LINE No.S ' 2300 2500 *,' 2900 * , 

'4600' AND ATTACH THE LABEL 'C' TO THE LINES, NUM“ 

BERED ' 2200 ‘ 2400 2800 ' AND '4500*. 
********************************-»<*-*^****'^**'^**’*^**'*^***** 


type 2a 
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APPENDIX-^ 

COMPUTER PROGRAMME No .6 


00100 C 
00200 C 
00300 C 
00400 C 
00500 C 
00600 C 
00700 C 
00800 C 
00900 C 
01000 C 
01100 C 
01200 C 
01300 
01400 
01500 
01600 

01700 9B 
01800 t: 

0 1 900 C 

02000 C 
02100 C 
02200 C 
02300 C 
02400 
02500 C 
02600 C 
02700 C 
02800 
02900 
03000 
03100 
03200 
03300 
03400 
03S00 
03600 
03/00 
03800 
03900 C 
04000 C 
04 1 00 C 
04200 
04300 
04400 
04500 
04600 
04700 
04800 

04900 100 

05000 C 

05100 C 

05200 C 

05300 

05400 

05500 

05600 

05700 

05800 


PROGRAMME FOR DRAWING THE GRAPH OF HEAT TRANSFER 
COEFF. VS. THETA FOR POLAR PROFILE TUBES, WHEN A VA- 
POUR IS CONDENSING OVER THE TUBES, USING THE INTERA- 
CTIVE GRAPHICS LIBRARY (IGL) SUBROUTINES. 

**'>i-*4f***********-»f************#************«-********** 

NOMENCLATURE: — 

X= ANGLE THETA, MEASURED IN RADIANS FROM THE UPPER 
GENARATRIX OF THE POLAR CURVE. 

Y«HEAT TRANSFER COEFFICIENT. 

P. 1,R.C.=PR0GRESSIVELY INCREASING RADIUS OF CURVATURE. 

*****#**#* •>t#'«-****#*->(-*****^***#*^»-*-i»-***->«-******'»*'***»^<-'*-** 

DIMENSION X(300) ,Y<300) ,AX(100) ,AY(10O> 

6022122 
DO 98 1=2,99 
X ( I) =X (1 ) ■^. 0314159* ( I-l ) 

CONTINUE 

********************************************** *««*#**♦*# 

THE FOLLOWING FILE READS THE HEAT TRANSFER COtTf . 

DATA FOR CIRCULAR PROFILE TUBES FROM AND T tit R PRO 
GRAMME, WHERE IT IS CALCULAFtl) AND STORED IN 1 HE 
SAME FILE. 

^^****4t********** ************ ♦^<***^***<»**^* ****** ^**^'^'**‘ 

READ(20,*) CY<n , 1 =^ 1 , 99) 

^*^^^<..j|.^(.4Ht.**4(.4(.^H*#********************'*'*^*****‘’*‘* 

INITIALISATION OF THE GRAPHICS TERMINAL. 

*->fr*-M-** ******** *********<***^ *********#****>»**«>* 

CALL GRSTRT (4010,5) 

CALL VWPORT (10.0, 11 0. 0 , 0. 0, 1 00* 0 > 

CALL W I NDOW ( - 1 . , 5 . 5 , - 1 000 • ♦ 14500 . ) 

XMIN=0. 

XMAX=5- 

YMIN=0. 

YMAX=14000. 

CALL M0VE(XM1N,YM1N) 

CALL DRAW (XMIN, YMAX) 

CALL MOVE (XMIN, YMIN) 

CALL DRAW (XMAX , YMIN) 

^^^i(.jn.+(.******************^ 

LABEILLING THE X-AXIS 

l(.^4^4t*****************-^* 

CALL TXICUR(8) 

DO 100 IX«=1,6,1 
AX<IX)»FL0AT<IX-1) 

CALL MOVEiAX (IX) ,0.0> 

CALL DRAW (AX (IX) ,-500.0) 

CALL MOVE (AXdX) ,-650.0) 

CALL RNUMBR(AX(IX) ,1 ,3*) 

CONTINUE 

**************^*****’*^** 

labelling the Y-AXIS 
**************'^^^'^**‘*^** 
call TXICURCl) 

DO 200 IY=1 , 14001 , 100® 

AY < lY) *=FLOAT ( lY- 1 ) 

CALL MOV£(0-0,AY<IY) ) 

CALL DRAW (-0. 1 , AY ( lY) ) 

CALL MOVE(-0.F,AY(1V)) 
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06000 200 
06100 C 
06200 C 
06300 C 
06400 
06500 
06600 
06700 
06800 

06900 10 

07000 C 
07100 C 
07200 C 
07300 C 
07400 C 
07500 C 
07600 
07700 C 
07800 C 
07900 C 
0B000 C 
08100 
08200 
08300 
08400 
08500 

08600 20 
08700 C 
08800 C 
08900 C 
09000 C 
09 1 00 C 
09200 C 
09300 
09400 C 
09500 C 
09600 C 
09700 C 
09800 
09900 
10000 
10100 
10200 

10300 30 

10400 

10500 


CALL RNUMBR(AY( IY> , 1 ,7) 

CONTINUE 

DRAWING THE GRAPH FOR CIRCULAR PROFILE TUBES. 

CALL MOVE ( X ( 1 ) , Y ( 1 > ) 

DO 10 J=2,99 
CALL SMOOTH 
CALL GRAIN (0.0) 

CALL DRAW(X(J) ,Y<J) ) 

CONTINUE 

************^t**************4t**»**<***->t***#->l-#-«-********** 

THE FOLLOWING FILE READS THE HEAT TRANSFER COEFF . 

DATA FROM ANOTHER PROGRAMME , WHERE IT IS CALCULATED 
AND STORED IN THE SAME FILE. (FOR P.I.R.C. PROFILE TUBES 
WITH BOTH GRAVITY AND SURFACE-TENSION EFFECTS. 

*■«•*■»»»■»#*#********#*##♦♦**#♦#«#*#■# #**♦«■**# 
READ ( 23 , * ) ( Y { I > , I 1 , 99 ) 

■^t-******')<-'H--H-#***-* *#*#**!-****# -I* « *«>** 

DRAWING THE GRAPH FOR P.I.R.C. PROF ILF LONS U>i F< I Nt5 
GRAVITY AND SURFACE-TENSION SlMUl lANFOUHlV. 

CALL MOVE (X ( 1) , Y ( l ) ) 

DO 20 K«2,99 
CALL SMOOTH 
CALL GRAIN(0.O) 

CALL DRAW(X (K) , Y (K) ) 

CONTINUE 

******»**->**'lt-«4t#* »**#« #»-»##*'IHH***#«iH**^»*'*«*»»*'»H»-**'««^'**^***'^ 

THE FOLLOWING FILE READS THE HEAT TRANBFER COEFF. DATA 
FROM ANOTHER PROGRAMME , WHERF IT IS CALCULATED AND STOP 
ED IN THE SAME FILE. (FOR P.I.R.C. PROF 1 1 E TUBf.S WUH 
GRAVITY EFFECT ONLY) 

**»##**«■*♦#* «**♦*»****►***** ••♦^ * « #« 

READ(24,*) (Y(I ) ,1-1 ,99) 

***'i<-*-)fr*->t-***##* #«»«»**« ♦*«'*********^***********‘*^'*** *“**** ** 

DRAWING THE GRAPH FOR P.I.R.C. PROFILE CONSIDERING 
GRAVITY ONLY. 

CALL MOVE (X (1) ,YC1) ) 

DO 30 M=2,99 
CALL SMOOTH 
CALL GRAIN (0.0) 

CALL DRAW (X (M) ,Y <M) ) 

CONTINUE 
CALL GRSTOP 
END 


. X 


.typ© u 
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A1.1-1m\1DIa -1o 


COoT 01" PKOiJUCTlON qk KLi'-iC'rR l CAI, i-'.N 
(SOURCE: PANKI THERMAL POVffiR STAl'lOH 


::kcy 

, KAN} -’UK 


/ 


Yr , of 
observation 


Year H^Tergy cost /unit [ P&ise_J 


1 . 

2 . 

3. 

4. 

5 . 

6 . 


1979-80 

32. '50 

1980-81 

35.67 

1 981 -82 

40.96 

1982-83 

46.24 

1983-84 

46.24 

1984-85 

57.50 

1985-86 

72.55 


7 . 
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APPENDIX -11 


00010 C 

00011 I 
00100 C 
00200 C: 

00300 c: 

00400 e: 
00500 C 
00600 
00700 C 
00B0W'i L 
00^^00 (' 

0 1 000 L, 

01100 C 
0 1 200 V 
0 1 300 C 
01400 (' 

01500 t; 
01600 t 
01700 « 

01B00 1 

0 1 «?00 
02000 ’>t 

02100 
02;.‘00 
0'2210 C 
02!>00 C 
02310 (’ 

02400 
02410 I 
02500 i ■ 
02510 I 
02600 
02 V 00 
02800 
027’00 I 
03000 
0 3 1 00! 

0*1200 
0:5*.00 
0 5401 ^ 5 

0.2700 i 
0 7600 
0*7 >*00 
0 7000 
0“',tV00 
04001 ^ 5 

04 100 
04 200 
04 701 * 

04 400 
04500 i 
04.700* < 

04V 00 10 

04B00 


COMPUTBR PRIK^RAMMI- No. 7 


PROGRAMME TO EVALUATE THE FUNCTIONAL RELATIONSHIP 
OF 'CEPU ' ,COSr OF PRODUCTION OF ENERGY WITH'TIME' 

BY EMPLOYING THE 'METHOD OF LEAST SQUARES'. 

********4(- *******♦*•«■ -tf **** 

D I MENS ION S ( 1 0 , 1 0 > , X ( 1 00 ) , & < 5 ) , A ( 5 ) , Y ( 1 00 > , C Y ( 1 00 > , PD ( 3 00 > 


NOMENCLATURE: 

M=NUMBER OF DATA POINTS 

N=.ORDER OF POLYNOMIAL TO BE FITTED THROUGH THE DA ! A POINTS 
X^INDEPENDENT VAR I ABLE ( T IME IN ' YEAR REPRESENTED BY 
THE No. OF OBSERVATION.) 

Y- DEPENDENT VAR I ABLE < ' CEPU ' ) 

CET-'U«COST OF F'RODUCTION Of- ENERGY. 

CY=T')EPENDENT VARIABLE CALCULATED FROM THE PDLYNOM I Al.-F 1 T . 
PDi-F’ERCEN T AGE DEV 3 AT I ON 
SUM=CUMUL,ATIVE 'PD ' . 


WRITE <49 ,31 ) 

F OKMA'I {///// /////// /'> 

H=:7 

THE^FOU.UriN^^^^ INPUT DATA TO THE ^ ^ ^ ^ ^ 

READ(24.H) <X <n ,Y( I) , 1 = 1 ,M» ^ ^ ^ ** »« 

EVALUAtIdN^ the ELEMENTS OF THF MATRIX 'S' ^ ^ 

SUMYs=0. 

SUMX Y- 0. 

SUM>:SY.-:0. 

SUMXI:Y-0. 

SUNX«0. 


BL)MXSR-*0. 
SUMXCB'"^0. 
SUM XI F7-0. 
SUMXFV=0. 


SUMXBX»0. 

DO 10 I«1,M 
SUMV=« S'UMN ( n 
SUMXY=SUMXY4-X (n*Y(I) 
SUMXSY-Sl !MXKY+ ( X H ) *¥7)^^ < I 1 
SUMXCY-SUMX(:Y'+ (X < I ) **3) ( I ) 

Sunx-STTru-t- X ( I > 
SUMXSR“SUMX£3R+X < I ) **2 
S I )M C B “• SI IH X t: H X ( I) * ^ 3 
sijMXFR-SUMXFW’*-X ( 1) +**4 
suMxr v~suMXFv-+ X ( n 
sutixsx=suMxsy.:^<n*«6 

tXJN I I NLIF: 

s o , n ~M 



12s 




8(1,2) -SlIMX 



S ( 1 , -8UMXS9 

iAV* t m 

( . 

5(1,4) -BOHXPH 

05;.' 01') 


S ( 2, 1 ) =-S8MX 

05 :.00 


5(:',2)--8llMX5K 

05400 


S(':',3)=-SLIMXCB 

05r..00 

t 

B ('2,4'--S8MXFR 

05600 


S ( 3., 1 )-SUMXSH 

05700 


S <; ‘■'.,2) -SUMXCH 

05B00 


S(3,3)=SUMXFR 

05900 

C: 

B ( 5 , 4 ) --SUMXf-' V 

06000 

c 

S <4 , 1 ) =SUhXCB 

06100 

i 

B(4 ,2)-SllhXKR 

06200 

c: 

S (4 ,3) -SUhXFV 


( 

B (4 ,4 )-'3UNXS>: 

064 00 


& (1 > =SUMV 

06500 


5(2) =^SUMXY 

06600 


& (3) ==SUMXSY 

06 7 00 

5 

H<4) •=-B8MxL V 

06900 


N="N4- 1 

06900 


WHl 1 E (21 ,->^) ( (5 ( 1 ,v1 ) , 1 ,W) , 1 1 . 

07000 


CLOSE (UN n -2 1 , DEV I CE== DSl , ) 

0 7 1 00 


Rt:9)D ( 2 1 , •► ) < ( B ( 1 , J ) , J ~ 1 , H ) , 1 - 1 , r 

07110 

( 


07 200 

t 

INVEKBlDf! OF- 1 HE MAIRIX 'B . 

07210 

V 


0/ :.00 


1 -• 1 

0/400 


N X “.N (• 1 

0“^1.00 


NY=-2-)*N 

0/600 


DU 80 J=NX,NV 

0 7 /00 


S ( 1 , J ) •'*• 1 . 

076100 


1=:H 1 

07900 

B0 

CUNT INUE 

08000 


L=^l 

0B100 


L “• 2 

08200 

1 10 

XM~S (L ,L ) 

08500 


DU 140 J=-1-,NY 

084 00 


G U ,..! ) (1. , J ) /XM 

08!:'0t'' 

140 

LOrU INUE 

08600 


DO 190 I'=K,N 

08700 


X !:-&< 1 ,l ) 

0BB00 


DO 191 J=L,Ny 

08901^ 


5(1 ,2)--=B( 1 ,2^-y(L ,J)*X3 

09000 

191 

CXJNTINUE 

09 1 00 

190 

UOl'l 1 INUF 

09 '2 00 


L/-'4 , -f 1 

09500 


1 — (. 4 1 

09400 


IF (L.-N) 110,1 10.230 

09500 

;*’/0 

L. = N 

09600 

27.f • 

1 /-'L-l 

09'H^0 


l)U 290 1 "• 1 ,l / 

0*7800 


84.-4 

09900 


’t t ■■ S (10 ■’ 

1 0000 

1 0 1 00 


DU 291 4*1 ,N> 

B ( 1 ,4 > V. ( 1 ,4 ) • *^ U ,4 / ♦►V 1 

10200 

2*^ 1 

Cl m } 1 Nt IF 

1 0’;.00 

290 

( UNI IWUl 

10400 

10500 

10600 

10/00 

5720 

1--4 -1 

U- U 1 > *' 

WRllh < <•.(! ..!>,.! >11-;. H.' ' , 

l;1 1 'St MINI 1 ' , 

, ii U-NX.N/**! 

10800 


READ(2.-,*) ( U.M » oJ MA , ^ 
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1 mw { 

1 1 vIVIU c 

1 li 01^' i 

1 1 
i 1 

1 1 40(.'> 

1 1 

1 1601*'' ■/ 1 
11700 /0 

1 1 800 l: 

1 1 W0 C 

1 0000 c 
12100 
1 :'200 
1 2300 
12400 I 
1 2500 C 
1 2600 C 
12700 C 
1 2800 
1 2900 
1 3000 
1 5 1 00 

a : '200 

13300 2R 
1 3400 
1 3500 
1 3600 

1 3 700 2 '/ 

1 '3800 
1 3900 

- t 


Mill lift ir.AlION ( 1( Tut iriVtI'Mti tUtll.!* lU WlHi ^U4U 1« M . 

DCl 70 1 = 1,N 

A ( n - 0. 

DC) 7 1 J•.= NX,^4V 

2 i < n = A < n ^ < I , .1 • f I' ' I f I ' 

CtlNI INHl- 
CUN ! INUl 

IHt ROW MA1RI7 n r, IVIR 1 HI- MH'IF.'. Of I Fil Fill! 2 C !1 . 


WR I T £ ( 4 B , ) •' A < n . I 1 , ( NN N M 

Cl l'S£ (UNI 1 - 4H,ni ^'ICf ■ ') 

RCAD (4H,* MA< 1 , I - I , <NV-N^ > 




CALCULAU(3N Of- THl: DRVIATIGN BHTWtbN 1 HE ACTUAl VAl.UH 
CAL cut A 1LU VALUE CH- IH& DEPENDLNl VARIABLE. 


AND 


SUM^0. 

DO 27 
CV (0 ) ^--0. 

DO 28 1=^1, N 

C; Y ( J ) =CY <CH'^(A(IH*(X(a)**(3"l>)) 
CONTINUE 

ET) ( J ) ” < a V ( J ) -• Y ( ;i ) ) / Y ( v'J > ) * 1 00 . 
SUM=SUM+F‘D(vD 

WR I I E ( 4V , # ) X ( J > .YU) , CV U1 ) , PI.) ( v1 ) , Bl IM 
CONI INUE 
S'i OP 
END 
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1 27 


Av'EfL'iGE HOUGfINESS OF COFiM.AjCLd, /~ 6_7 


Material Absolute rougliness of pipe 

( in mm.) 


Riveted steel 


0.9 -9.0 

Concrete 


0.3 - 3.0 

\/ood stave 


0.18- 0.9 

Cast Iron 


0.26 

Galvanized Iron 


0.15 

Asphalted Cast Iron 


0.12 

Commercial Steel or 

V»'r ought 
Iron 

0.046 

Drawn Tubing 


0.0015 


Glass 


Smooth 


APPEND 1X^1 3 


m 1 c 

00200 C 
00300 C 
004190 C 
00500 C 
00600 C 
00700 C 

00800 r. 

00900 i; 
01000 C 
01100 f: 
0 1 200 C 
01300 (■ 

01400 C 
01500 r: 

01600 r 
01/00 {• 
0 1 800 (/. 
01900 I 
02000 C 
02100 L 
02200 C 
02300 t: 
02400 C 
02500 C 
02600 C 
02700 c: 

02800 C 
02900 c: 
03000 i: 
03100 L' 
03200 C 
03' •.00 (• 

03400 r: 
03is:.00 i 
03600 C 
03700 L 
03800 C 
03900 t.. 

04000 r 
04100 I 
04200 i: 
04300 ( 

04 400 (' 

04 5.00 l 

04600 r 
04 /00 t 
04800 (. 

04V00 I 
05000 C 
0!»100 I 
05200 L 


COMPUTER PR(3BRAMME No. 3 


•» ***-»(*-».# If 

f-‘RGURAMME F OR GH'l A 1 N 1 NG 1 HE Df-U 1 Mt m C (INDE rsJS'.EP 
TUBE DIMENSIONS ON THE BASIS Of- MINIMUM HU Ai 
C 0 S 1 F UP C I RCUl . AK 'I I JHE PR DF I L t . 

NOMtNCl AlUKT 


AI^'-^^AREA FACTOR 

MUL=DYNAMIC VISCOSITY OF 1 T fXH O C()NI)CNS/UF 
FI -^THERMAL CONDUC ( 1 V IT V OF I , J UU I D Of »N1)F F 

L..r»*LENGTH OF TUBES 

C:F:F’U-~C0S*1 of F'RODUf'TinM df^ fnfri^y 
HT CF ==HE A 1 T RANSF F F‘ COEF F 1 C 1 EN 1 F AL U IR 
QOR^ORIGINAL HFA'l TRANSFER RA I E IN THE rUNOANSFR 

etap^overall efficiency of the pump 

DO^OUTER DIAMETER OF TUBES 
Di=INNER DIAMETER OF TUBES 
EPSN-ABSDL LITE ROUGHNESS DF TUBES 
INTRT=RArE OF INfEREST' 

RHOL-MASS- DENSn Y OF THE LIQUID CONDENSAU- 
LMTD=LOG MEAN TEMPERATURE DIFFERENCE 
PR=.PRANLV1L No. 

HS=STEAM~S1DE HEAT TRANSFER CGEFFICIENV 
RHOM=MASS-DENSITY OF MATER I AL , USED I!4 THE TUBES 
CM=COST OF MATERIAL. 

VFR~ VOLUME FLOW RATE OF COGl.ING WATER IN T HF CDNDF NSER 
PER= PERIMETER OF TUBES 

acn~arfa of cross section 

RR=RELATTV£ ROUGHNESS OF TUBES 
VEl -VELOC ITY OT L'.OQLING WATER 
HMD-HYDRAUl. IC MEAN DIAMETER 
REN-REYNOl D'S No. 

FF =FR I C T I ON FACT OR 

HW-WATER*'B1DE heat transfer COEFF Jf/IENT 
THRS-THERMAL RESISTANCE OF STEAM SIDE 
THRW^OHETRMAL REFUSTANCF OF WATER SIDE 

thrm*-.mis(:el,i.emeou5 thermal resistance 

I H t- O'v'Efy Al L HE AT T RANSF E R f 5DE FF- U' I E NT , BABE D 
ON OUTER SURFACt AREA 
A(:i=* OUTER BLIPI'ACE AREA 

DP..^PRESSURE DROP OF COULIN8 WATER IN CUNDENSR lUBES 

DYNF’^UYNAMJ C PRESSURE OF COOT ING WAI F.R 

PWRT = THEOR I T I CFU- F*l IMP I NU PUWf: R 

pi*jpp=;,AC T I ir.l PUMf ' 1 f-iU PI iWF R 

EWCV-ENERI.. r' l.UN'.UMPT ION PFR 

C.L N-COST 14 FN!T-u.' 

CFNT-TOTAL COST <4 fcNLRi,f 

F'W£' NO-pRt f:i N I W if* U l Of F NMn U 11 C. I 

VOL M-*VOLUMT Of Mm r I K I Ai . » T'>T T' 
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10/00 
10900 
10900 
1 1000 
11100 
1 1 200 
1 1300 
1 1400 
1 1 f.00 
1 1 600 
1 1 700 
1 1800 
1 1 900 
12000 
12100 
12200 
1 2300 
12400 
1 2500 
12600 
12700 
12800 
12900 
13000 
13100 
1 3200 
13 700 
1 3400 
1 3500 
1 3600 
13700 
1 3800 
1 3900 
14000 
14100 
3 4 200 
14300 
14400 
14500 
i 4600 
1 4 700 
14800 
1 4900 
ir.000 
15100 
15200 
15300 
15400 
1 5500 
1 5600 
15700 
15800 
15900 
16000 
16100 
3 6200 
16300 
16400 
1 6500 
1 6600 
16/00 
1 6B00 


9 


l: 

C 

(; 

7 

10 


i: 

c: 

c 

C' 


19 


29 


C 

c 

r 

t, hHb 

c; 

i, 

c 

C: 


UU=1 . / SUM 7 HR 
A0= (QOR) / <UO*LHT D) 

1. 1 UAO) / (PI ^ DON ( n ) 

DF ~(FF)*(l T)*( VF t ‘ I }vj ^is 11 > i' ^ «( » i ; i , , 

!)YNF-- (0.5)# (FclUH ^ # < VI i 

DP=DYNP+DP 

PWF;*] =. ( DP) ^ ( VFK ) 

PWRP=: (PWRT) / (F I Af ) 

E NL Y - ( PWRF' ) •<» ( 24 . 0 > » ( ' ( 2 . . 0 . / ( j t’i00 0 ■ 

CENT ==0.0 
DU 9 J=S,32 
X=FLOAT (J) 

CEN«t:NUY*CEP( 1 ( V I / ! 00 . 0 
PWENC=:CEN / ( ( 1 . 0 * I H 1 P 1 > # »« ( ;| fn 1 
CENT wCENT+PWENt: 

CONTINUE 

'70tM=:(3. f .Dlifi. f . i ((UU*# .»*P ti 

MTM« (VOLM) * (RHUM' 

LMT = ( M T M ) ^< ( (,:M ) 

CT (I )*CENT4-CMr 

THE F0L1.0WIN(3 MIF STUPFl. IMP IS Of I Hi ADHVE ANfU 

WRITE ( 47 , 7 T DON ( I ) . f f , FfJi, V , CPU f ,rM T ,0 ! ( [ ) ,M I M , t 1 
FORMAT ( /5X ,F5. 3 , 3>: , f 1 0. R , , f- 1 5 . 8 , 3X , f 31 . 2 , i B . 2 , 

1 3X , F 1 1 . 2 , 3X , F / . , 7X , F 6 . 2 ) 

CDNT INUl 

STOP 

END 


SUPKGIJ T 1 Nf: fur dr ( a I N I NO F R if: T I ON F AC TOP F RDM 
THE IMPl. 3C:ll Ri'IATlON HY ' CDl E BROUi. ' i.lSH'-JC: NF.*- 
W T ON-RAPHSON T E CHN I DUE . 

SUBROUTINE FR I XCJN <RR , RFN , FF ) 

X=0. 005 

DF X 1 =2 . 5 1 / ( REN-«-BOf< f ( X ) ) 

DFX2--RR/3. 7 

DFX3«DFX1+DFX2 

F X 1 -- ( -1-0. 86859 ) * ( Al 08 ( Di X.3 > ) 

FX2=1 . /SORT (X ) 

F X=FX3 +FX2 

DFXA~ ( 1 . /DFX.3) « (0. B6859) 

DF XB~ ( - 1 . 255/f<F N ) / o' # > 1 . 5) 

OFXL‘=^-0.5) /(X^<#1.5) 

DFX- (DFXA^DF XD) -4 ( Di XD) 

XN«<X)-<FX/DFX) 

I F ( ( ABS ( X'N- X ) ) . I F . 0 . 00000 3 ) 00 1 O ?9 


X=:XN 

on (0 19 
FT = XN 
KT- rUF<N 
END 






NOTE: — 

(0 USE THE- SAME PROGRAMME fOlP PUl AR PROFilE*-. 
OTHER THAN CIRCUL Af( , 1 HE VALUES OF ' AF E^ND 'HTCF' 

1 f I HE- CHE^NC:!' D SI. « 1 ( AHl Y , AS D I SCi IBSF D I M H il 

■CHAPTER. 2’ OF THE THEG I 'S . HOWEVER , THE CORRE SFTTND I N( > 
VmIUES of fHH r.HHVE PARAMETERS ARE flVEP AS 'UNIfV; 
FOR CIRCULAR PI OFILE: TUHEJ: 


. t ypr 



REFERENCES 


G.G. Shklover, V.p. Semenov, V.V. Pryakhin and 
A,M, Usachev, Heat Transfer with Steam Condensing On 
a. Horizontal Tube with a Profile Of Variable Curvature *1 
“Hiermal Engineering, vol, 32, No. 3, pp. 125-127, 1985. 

S.W, Churchill, Laminar Film Condensation** , International 
journal of Heat and Mass Transfer, vol. 29, No. 8, 
pp. 1219-1226, 1986. 

Frank Kreith and Mark S.Bohn, ^‘^Principles Of Heat Transfer"**, 
4th edition. Harper International Edition, New -York, 
pp. 542-552, pp. 646-652, 1986, 

W.F. Stoecker and J.W. Jones, Refrigeration and Air 
Conditioning**, 2nd edition, McGraw-Hill, Singapore, 
pp. 245-252, pp. 416-417, 1982. 

W.M. Rosenhow and J.P. Hartnett, eds., ^^Handbook of Heat 
Transfer**, McGraw-Hill, New York, pp. 12.1- 12.19, 1973. 

F.M. White, Fluid Mechanics**, 2nd edition, McGraw-Hill, 

New York, pp. 308-316, 1986, 

ft 3) 

O.G. Titejens, Fundamentals of Hydro and Aero -Mechanics , 
Dover Publications, New York, pp. 60-64, 1934. 

J.R. Welty, C.E. Wicks and R.E. Wilson, ^Vtmdamentals 
of Heat, Mass and Momentum Transfer** , 3rd edition, 

John Wiley and Sons, New York, pp. 390-397, 1984. 



132 


9. 


10 . 


11 . 

12 . 

13. 


14. 


15. 


16. 


V. Rajaraman, Computer Oriented Numerical Methods” , 

2nd edition, Prentice Hall of India, New Delhi, 
pp. 132-144, 1985. 

S.A. Hovanessian and L.A. Pipes, Vgital Computer 
Methods in Engineering”, McGraw-Hill, New York, pp.1-34, 
pp. 145-155, 1969. 

V. Rajaraman, Principles of Computer Programming ^ , 

2nd edition, Ek*entice Hall of India, New Delhi, 1985. 

PLOT 10, 4010C01, TEK User’s Manual , Interactive 
Graphics Library, 1982. 

Anthony J . Pettofrezzo, Introductory Numerical Analysis , 
D.C. Heath and Company, Boston, pp. 51-58, 1967, 

S.Kakac and T.N. Veziroglu, eds., Two -phase -flow and 
Heat Transfer , Hemisphere Publishing Corporation, 
Washington, pp. 817-842, 1977, 

f € 

Xue-Jun Chen and T. Ne^jat Veziroglu, eds.. Two-phase 
flow and Heat Transfer , Hemisphere Publishing Corporation, 
Washington, pp. 283-307, pp. 355-357, 1984. 

ts 

D. Buttei»worth and G.F. Hewitt, eds,. Two-phase flow and 
Heat Transfer , Oxford University Press, Oxford, 
pp. 459-462, 1974. 

** h 

W. Nusselt, Die Oberflachenkondensation des Wasserdampfes , 

Z. Ver. dt. Ing. 60 , pp. 541-546, pp. 569-576, 1916. 


17. 



133 


18. C.p, Gupta and. Rajendra Prakash, ^^Engineering Heat- 
'Transfer , 4th edition, Nem Chand and Bros,, Roorkee, 
pp. 346-354, 1983. 

19. Frank P. Incropera and David P. Dewitt, ^‘Fundamentals 
of Heat and Mass Transfer , 2nd edition, John Wiley end 
Sons, New York, pp, 480-489, 1985. 

20. P.J, Marto, D.J, Looney, J.W, Rose and A.S. Wanniarachchi, 
(€ 

Evaluation of organic coatings for the promotion of 
dropwise condensation of steam V International journal of 
Heat and Mass Transfer, vol. 29, No. 8, pp. 1109-1117, 
1986. 

21. Thomas F. Irvine, Jr. and James P. Hartnett, ads.. 

Advances in Heat Transfer , vol, 9, Academic Press, 

New York, pp. 181-268, 1977. 



MB - 1^07 -rvi-fAT-^ Con 



